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ABSTRACT 

We present a new study of stellar mass in a sample of ~ 70 submillimeter-selected galaxies (SMGs) with 
accurate spectroscopic redshifts. We fit combinations of stellar population synthesis models and power laws to 
the galaxies' observed-frame optical through mid-IR spectral energy distributions to separate stellar emission 
from non-stellar near-IR continuum. The availability of spectroscopic redshifts significantly enhances our 
ability to determine unambiguously not only the mass and luminosity of SMGs, but also the presence and 
contribution of non-stellar emission to their spectral energy distributions. By separating the stellar emission 
from the non-stellar near-IR continuum, we find that ~ 50% of our sample have non-stellar contributions of less 
than 10% in rest-frame //-band and ~ 10% of our sample have non-stellar contributions greater than 50%. We 
find that the ZT-band luminosity of the non-stellar continuum emission is correlated with hard X-ray luminosity, 
indicating an AGN origin of the emission. Upon subtracting this AGN-contributed continuum component from 
all of the galaxies in our sample, we determine a lower median stellar mass for SMGs than previous studies, 
~7x 10 10 M Q . We use constraints of the starburst time-scale from molecular gas studies to estimate the amount 
of fading our sample would undergo if they passively evolve after the starburst terminates. The results suggest 
that typical SMGs, while among the most massive galaxies at z ~ 2, are likely to produce descendants of similar 
mass and luminosity to L* galaxies in the local universe. 

Subject headings: galaxies: high redshift — galaxies: formation — galaxies: evolution — infrared: galaxies 



1. INTRODUCTION 

Soon after their discovery, the population of high-z, far- 
infrared (far-IR) lumi nous galaxies revealed in deep submil- 
limeter surveys (e.g.. ISmail et al.|[T997t iHughes et al.l fl998: 
iBarger et ai]|199&lEales et alj|1999l7 were proposed to repre- 
sent the formation of (the most) mass i ve, metal-rich spheroid s 
(lLillv etal.1 [19991 ISmail et alJ 12001 LSwinbank et all [2006). 
Submillimeter-selected galaxies (SMGs) have thus tantalized 
observers as a potential opportunity to witness the forma- 
tion of massive galaxies directly. The far-IR luminosities 
(> 10 12 L W ; e . g.. IChapman et all 120051 iKovacs et all [2006; 
iMagnelli et alj 120101: IChapman et alj 1201 Oh of these dusty , 
gas-rich galaxies (e.g.. IFrayer et al.lll998t IGreve et al.ll2005l) 
imply star formation rates in excess of lO 3 M yr _1 , suf- 
ficient to form the stellar mass of an L* galaxy in ~ 
lOOMyr, potentially explaining the seemingl y uniformly- 
old st e llar populations of loca l ellipticals (e.g., iBower et alj 
[1991 iKauffmann etail 12001 and ma ssive g alaxy popu- 
lations at z ~ 1-3 (e.g.. iFranx et alj 120031: ISteidel et alj 
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l200l lGlazebrooket~aT] Hooi iDaddi et all l200l. More- 
over, the observe d clustering of SMGs (Blain etal. 2004b; 
WeiB "eTaTJ 120091) is similar to that of massive galaxies in 
the local universe, and suggests that SMGs reside within 
~ 1O 13 M dark halos. Yet, the formation of such mas- 
sive galaxies on such rapid time-scales is difficult to ac- 
count for using theoretical models of galaxy formation 
(e.g.. iGuiderdoni et al.lll998t iDevriendt & Guiderdonl 120001; 
Bau gh et al.ll2005USwinbank et alj|2008h . 

Determining the masses of high-z SMGs, especially the 
stellar and gas masses, and their relationship to other popu- 
lations of high-z galaxies is thus crucial to understanding the 
link between SMGs and massive galaxies in the local Uni- 
verse and places much-needed constraints on the z = de- 
scendants of SMGs. Prior to the advent of the Spitzer Space 
Telescope, however, constraints on the stellar mass were dif- 
ficult to obtain for SMGs, especially for individual galaxies, 
due to significant obscuration at optical and near-IR wave- 
lengths (rest-frame UV and optical; see ISmail et aT] 120041 ; 
ISwinbank et~aT]|2004l; IPope et al.ll2005l) . The rest-frame near- 
infrared (near-IR), which Spitzer gives us the power to probe 
for high-z galaxies, is less affected by reddening than optical- 
band data, and thus is beneficial when looking at SMGs which 
suffer from extinction in the rest-frame UV and optical bands. 
In addition, because the near-IR luminosity is much less de- 
pendent on past star formation history than, for example, th e 
optical-band luminosity (e.g., IKauffmann & Charloll 1 1998b . 
the rest-frame near-IR continuum emission permits the best 
determination of total stellar mass in bot h young and ol d stars . 

Using data f rom Spitz e r-IRA C, lBorvsetal.1 d 2005), 
lHainlinel d2008h . iDve et all d2008l) . and Michaiows ki et alj 
(2010) have all calculated stellar masses for individual SMGs 
by fitting the observed-frame optical through mid-IR spectral 
energy distributions (SEDs), obtaining typical stellar masses 
rangin g fro m 10 10 8 - 10 u 8 M ? ; lHainlinel J20081) . IDve et alj 
(2008), and lMichalowskietalJ ( 1201 Oh also find evidence of 
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significant stellar mass assembly prior to the SMG phase. 
However, none of the previous stellar mass studies addressed 
the contribution of near-IR continuu m ex cess to th e SEDs of 
SMGs, noted by IBorys et all (120051) and iHainlind (120081) in 
numerous SMGs in their samples, which may represent very 
hot dust heated by AGN. Contamination of the stellar masse s 
by AGN could especially be a problem for B orvs et ail d2005l) . 
who used a small, X-ray-detected SMG sample which was 
thus possibly biased towards SMGs with stronger AGN con- 
tributions (or lower X-ray obscuration). While X-ray observa- 
tions and mid-IR spectra of SMGs indicate that the majority of 
the population are not dominate d by an active nucleus (AGN) , 
many SMGs host AGN activity dAlexander et alj|2003l 120051: 
iPope et al.l 120081; iMenendez-Delmestre et all 120091) . which 
may not be negligible for the purpose of stellar mass studies. 

We have published a Spitzer-IRAC and MIPS survey of a 
sampl e of ~ 70 radio-dete cted SMGs with spectroscopic red- 
shifts (Hai nline et al.ll2009l hereafter Paper I), a factor of 5 
larger than the sample considered by IBorys et al.l (120051) . In 
this paper we combine our IRAC data set from Paper I with 
optical and near-IR photometry from the literature to analyze 
the rest-frame UV to near-IR SEDs of the sample to determine 
stellar masses for these SMGs, while accounting for a con- 
tribution to the near-IR luminosity from hot, possibly AGN- 
heated, dust, expanding the analysis of Hainline (2008). With 
our large sample size, spectroscopic redshiftjTand treatment 
of non-stellar emission, we can make estimates of SMG stel- 
lar mass which better represent the stell ar component than 
tho se estimates curren t ly ava ilable from IBorys et aTl (120051) 
and iMichalows ki et al.l (120101) . Using our new estimates, we 
can evaluate whether SMGs contain enough mass to represent 
the formation of the most massive galaxies and make predic- 
tions for SMG descendants. 

The structure of the paper is as follows: we begin in $2l 
describing our SMG sample and data. In ^3]we discuss the 
methods we use to determine stellar masses for our sample. 
In ^4] we present our results from the fits of the stellar popu- 
lation models to the observed SEDs. In $5] we conduct inde- 
pendent checks on the stellar masses from dynamical masses 
obtained through observations of CO rotational transitions, 
with the caveat that the dynamical masses are also subject 
to systematic uncertainties arising from assumptions that the 
gas traces the total mass and the galaxies are in virial equi- 
librium, then compare the stellar masses of SMGs to those of 
other high-z galaxy populations, and conclude by discussing 
the evolution of SMGs to the present day. Throughout our 
analysis and discussion, we assume fijyi = 0.27, Oa = 0.73, 
and H = 7 1 km s" 1 Mpc" 1 dHinshaw et al.ll2009l) . 

2. OPTICAL, NEAR-IR, AND MID-IR DATA 

In Paper I, we present Spitzer-IR AC measurements for 
the spectroscopic SMG sample from [Chapman et al.1 (I2005L 
C05 hereafter). Their final spectroscopic sample contains 
73 SCUBA galaxies (median S^so = 5.7±3.0mJy), each of 
which was required to be detected in very deep (la noise 
~ 5 — 11 /iJy), high-resolution (~ 1 ") Very Large Array ( VLA) 
maps at 1 .4 GHz to obtain a position accurate enough for op- 
tical spectroscopy. The redshift distribution of the C05 sam- 
ple, when corrected for spectroscopic incompleteness, is well- 

1 1 Note that using optical through near-IR SED fitting to obtain both pho- 
tometric redshifts and stellar mass for SMGs, as done in Dye et al. ( 2008), is 
problematic in that both the redshift and mass are sensitive to the choice of 
spectral template used. 



characterized by a Gaussian with median redshift of z = 2.2 
and er z = 0.8. None of the z > 1 SMGs in the final sam- 
ple analyzed in C05 are suspected to be strongly lensed. 
We refer the reader to C05 and Paper I for further details 
on the sample and Spitzer observations. We will exclude 
from our analysis 2 SMGs from the IRAC-observed sub- 
sample whose radio counterparts are in doubt (and which 
were not detected by Spitzer): SMM J 1 3 1 23 1 .07+424609.0 
and SMM J221725.97+001238.9. We also exclude from our 
analysis the galaxies which are suspected to be low-z lenses 
of background submillimeter sources in C05, as in C05 and 
Paper I. The final sample which we analyze in this paper 
thus contains 67 radio-detected SMGs, and includ es the entire 
sample of 13 X-ray detected SMGs analyzed by Borys et al. 
(2005). We note that a bias towards emission-line galaxies is 
introduced into our sample since the SMGs are required to be 
spectroscopically-identified: the redshift of an emission-line 
galaxy is more easily identified than that of a galaxy without 
emission lines when the continuum flux is faint and has a low 
signal-to-noise ratio. While one obvious consequence of this 
bias is that our sample is incomplete in the "redshift desert" 
between z ~ 1 .2 and z ~ 1.8, the bias may also result in pref- 
erential selection of galaxies which contain less-obscured star 
formation or AGN activity. In addition, since we require a 
radio detection of the galaxies in our sample in order to ob- 
tain the spectroscopic redshift, we may be selecting against 
SMGs with cold dust or higher redshift (see, e.g., Ivison et 
al. 2002; Eales et al. 2003; C05). As a result of these biases, 
the conclusions we draw must be restricted to apply only to 
spectroscopically identified, radio-detected SMGs. However, 
a wide range of luminosity and redshift is still probed within 
these s election limits (see, e.g.. lIvison et al.l l2002; Blai n et al] 
l2004ah . 

Although we postpone a discussion of the uncertainties as- 
sociated with our modeling procedure to later sections, we 
note that, in principle, to obtain the best constraints on the 
stellar population ages and the extinction of starlight from fit- 
ting stellar population models to the SEDs of high-z galaxies, 
we must use photometry which spans the 4000 A break and 
the rest-frame UV t hroug h optical bands, respectively (e.g., 
Shapl evet alJl200lT 120051) . Thus, for our present study we 
combine observed-frame optical and near-IR measurements 
from a variety of literature sources with our IRAC measure- 
ments presented in Paper I to obtain the most wavelength cov- 
erage possible. The available observed-frame optical photo- 
metric b ands vary for each of the seven sky fields in our sam- 
ple; see jj3.3l for a discussion of the impact on our final results 
of the varying number of photometric constraints within the 
galaxy sample. In Table Q] we summarize for each field the 
photometric bands available and the typical number of opti- 
cal and near-IR photometric constraints per galaxy. For the 
CFRS -03h field, we use here UBV data from Clements et al. 
(2004) and R data from C05. For the Lockman East field, 
we use the B data when available from C05 and R data from 
llvison etalJ d2005b . In the HDF/GOODS-N field (hereafter 
simply the GO ODS-N field), we use the optical photome- 
try catalogue of Cap ak et alJ (|2004), which contains measure- 
ments in U, B, V, R, I, and z'. The SMGs in the SSA-13 
field have B an d R measurements in C 05, plus z band mea- 
surements from Fo malont et al.l ([2006) of comparable depth. 
In the CFRS-14h and SSA-22 fields we use the B and R mea- 
surements from C 05. For the ELAIS -N2 field, we take our 
optical data from llvison et"aTI (120021) . in B, V, and R. The 
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available near-IR data for our sample are more homogeneous, 
because we obtain all of th e near-IR data (/, J, K) from the 
study of S mail et ail (120041) for all fields with the exception of 
the G OODS-N field For SMGs in GOODS-N, we use the 7 
data oflCapak et al.l (120041) . but the J and K data o flSmail etail 
(2004). We convert all the optical data obtained from the lit- 
erature to total magnitudes using aperture corrections, gener- 
ally supplied by the corresponding authors of each data set, to 
eliminate the effects of differing seeing conditions and pho- 
tometry aperture size between different data sets. 

We note that 4 out of the 67 SMGs from our IRAC-observed 
sub-sample of C05 SMGs have less than 3 detections across 
the observed-frame optical-mid-IR wavelength range. While 
it is difficult to place meaningful constraints on their stellar 
SED to derive ages and extinction with so many upper limits, 
we can still place upper limits on their luminosities in a given 
band, and use mass-to-light ratios derived for other SMGs to 
place upper limits on their stellar mass. Thus, we retain these 
galaxies in our sample despite the small number of detections. 

3. CALCULATING THE STELLAR MASS OF SMGS 

Our approach to deter mining the stellar masses of SMGs 
is derived from t hat o f iBorvs et al . (2005), because, as we 
discuss below in §13.11 we find that the method typically ap- 
plied to obtain the stellar mass of intermediate- and high- 
redshift galaxies (i.e., fitting a stellar population model to each 
galaxy's SED to determine its stellar population age, visual 
extinction, and appropriate mass-to-light ratio) is unable to 
adequately constrain the stellar mass of individual SMGs in 
our sample. Rather than assume that the stellar population pa- 
rameters and masses derived for individual SMGs from SED 
fitting are correct in an absolute sense, we derive constraints 
on the stellar population parameters and stellar masses from 
averages over the SMG sample ( $3. 2\ . In doing so, we can 
only constrain the stellar mass of SMGs in an average sense 
(the mean and median of the sample). Although we will still 
present estimates of stellar mass for the individual galaxies for 
those who wish to use them despite the large uncertainties, we 
will only discuss the average properties of the sample. 

In the discussion that follows, and to obtain the average 
stellar population parameters of SMGs, we have used version 
1 1.0 of the HYPER- Z p hotometric redshift software package 
dBolzonella et al.l [20001) to fit evolutionary population syn- 
thesis models to the optical-mid-IR SEDs of the individual 
galaxies in our sample. HYPER-Z employs % 2 minimization 
to derive the best-fit age/time since the start of star formation, 
optical extinction A v , and normalization for a given stellar 
population model. We restricted the code to fit the SEDs at 
the spectroscopic redshift of each SMG, and required that the 
best-fit ages be less than the age of the Universe at the redshift 
of each S MG. We have fitted t he p opulation synthesis mod- 
els of both Bru zual & Chariot! (120031 using the Pad ova 1994 
stellar tracks; hereafter BC03) and Maraston (2005, hereafter 
M05), although in §|4] we choose to report results for only 
one set of models. We u tilized the solar metall icity models 
based on the findings of ISwinbank et al.l (120041) that SMGs 
have metallicities consistent with those of UV-selected galax- 
ies of similar redshift, which for the most massive exa mples 
have approxima tely solar metallic ities (Erb et al. 2006aJ), and 
the findings of iTecza et all (120041) . The BC03 and M05 mod- 
els are made available with two different initial mas s func- 
tions ( IMFs): the BC03 models hav e a choice of the lSalpeter] 
(fl955l) IMF or the IChabrieri (2003) IM F, while for the M05 
models we must choose between the ISalpeterl (119551) IMF 



and the Kroupal d2001l) IMF. We preferred to not use the 
IMF of Salpeterj (119551) since observations of local galaxies 
suggest that this IMF over-predicts the num bers of low-mass 
stars and thus the ma ss-to-light ratio (e.g.. iBlain et aljfl999t 
Bell & deJong 2001), casting suspicion on its use in high-z 
galaxies. Instead, we use the[Chabrier (2003) IMF for the 
BC03 models with a lower mass cutoff of 0. 1 M w and upper 
mass cutoff of 100M Q , and the iKroupal (120011) IMF for the 
M05 models with the same low-mass cutoff. Although the use 
of two different IMFs may be somewhat confusing, the two 
different IMFs predict similar spectral properties and mass- 
to-light ratios at a given age (BC03), and lTacconi et al.l ( 120081) 
find that the two IMFs produce similar results to within 15%. 
To account for dust extinction in our SED fitting, we assumed 
a simple dust screen model, although it is likely that SMGs 
contain a variety of clouds of different obscurati on attenuat- 
ing th e total observed starlight. We used the Calze tti et al.l 
(2000) extinction law for starburst galaxies, allowing the ex- 
tinction to vary within a range of Ay = 0-5 mag. In addition 
to extinction intrinsic to each galaxy, we also corrected for the 
modest reddening along the line of sigh t to each SMG due t o 
the Milky Way using the dust maps of iSchlegel et al.l (119981) . 
The line-of-sight reddening estimate for each sky field is listed 
inTableQ] 

3.1. Complications In Applying Standard SED Modeling 
Techniques to Estimate M^ for SMGs 

3.1.1. Constraining Star Formation History of SMGs 

For our stellar mass study, we had originally intended to al- 
low HYPER-Z to choose the best-fitting star formation history 
for each SMG from a set comprised of an instantaneous star- 
burst star formation history (which is equivalent to a simple 
stellar population or SSP), a continuous star formation his- 
tory, and exponential star formation histories of the form 

SFR(f)ocexp(-f/T), (1) 

where r = 0.1,0.25,0.5,1,2, and 5Gyr, and thus obtain an 
unambiguous stellar population age, Ay, and stellar mass-to- 
light ratio which could be used to calculate the stellar mass of 
each galaxy. However, careful inspection of our tests utiliz- 
ing the various star formation histories reveals that for many 
SMGs in our sample, the fits could not indicate a unique best- 
fit star formation history. We find that the values of xl cor- 
responding to the best fit Ay, age, and SED normalization are 
extremely similar between even just this restricted set of star 
formation histories, as are the appearances of the fits. Ad- 
mittedly, for some SMGs in the sample with relatively few 
photometric constraints or anomalous observed-frame optical 
photometry, such a result is not unexpected; however, we en- 
counter the same ambiguity among galaxies with the most 
complete and accurate photometric data. We illustrate the 
problem in Figure [T] where we show the observed SED of 
a typical SMG at z ~ 2.5 with excellent photometry, fit by a 
range of different simple star formation histories, using mod- 
els of both BC03 and M05. Despite the broad range of pho- 
tometric data in Figure Q] {UBVRIz' JK+lRhC), nearly all of 
the star formation histories can produce fits of similar like- 
lihood (e.g., 4 of the 6 histories shown produce xt values 
within ±0.01 of each other for the BC03 models, and all of 
the xl values correspond to a probability of 40-70% of ex- 
ceeding x 2 for the number of degrees of freedom), with ages 
varying from ~ 10 Myr to ~ 2 Gyr and extinction in the range 
A v = 1.5 -2.4 mag. 
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Figure 1. Fits of several different stellar population models of different star formation history to the observed SED of SMM J 123707.21+621408.1 (round points) 
at z = 2.484. The best-fit BC03 model spectrum is over-plotted as the solid black line, while the best-fit M05 model is represented by the dotted (green/gray) line. 
The different fits are nearly equally acceptable in a \ 2 sense. Note that this particular SMG does not display a near-IR continuum excess. 



Thus, we conclude that we are unable to place use- 
ful constraints on the star formation histories of individual 
SMGs through fitting such simple stellar population mod - 
els, asJsFaplev et alT d200Tl 120051) . iPapovich et all J200l . 
and lErb et alj d2006bl) also find for UV-selected high-redshift 
galaxies with similarly broad wavelength coverage. Since 
the best-fit ages are strongly dependent on the assumed star- 
formation history (the instantaneous burst models tend to 
produce the youngest ages, while the continuous star for- 
mation models tend to produce the oldest ages, see, e.g., 
Maras ton et alj |2010), we are also unable to effectively con- 
strain the ages of individual galaxies and thus the appropri- 
ate mass-to-light ratio of the stellar populations in individual 
SMGs without other information on the likely star formation 
histories, introducing extra systematic uncertainties into the 
stellar mass calculations, on top of any resulting from the stel- 



lar models themselves. 

Some of our inability to constrain the star formation his- 
tory and the stellar population ages stems from the avail- 
able photometric data for SMGs. A small subset of SMGs 
in our sample (10 out of 67) possess seemingly anomalous 
observed-frame optical photometric data; for these galaxies, 
the universally-high xi values indicate that none of the stel- 
lar population models adequately describe the SEEQ Yet, 
even when the wavelength coverage and data quality are ex- 
cellent (e.g., the GOODS-N field), we still may have insuf- 
ficient information. The determination of age from stellar 

12 The causes of the anomalous photometry are unknown, as no notes are 
available from the literature data sources for these galaxies. In a few cases, 
the rest-frame wavelengths of the discrepant data suggest a contribution from 
an emission line; however, in the remaining cases the optical counterpart 
could have been mis-identified or the photometric error is larger than stated. 
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population models relies heavily on the shape of the Balmer 
break or 4000 A brea k in the SED (see, e.g. lKriek et"aT1l2008t 
IWhitaker etdl ^QlQlf^l Unfortunately, for only up to a third 
of SMGs in our sample do we see a suggestion of a break 
in the continuum near 4000 A, and we usually do not have 
enough photometric data points for our SMGs in the criti- 
cal area of the Balmer/4000A break to clearly delineate the 
location and the shape of the spectral break. In addition, 
the well-known degeneracy between age an d extinction (e.g., 
ISawiclri& Yeelll998t IShapleyet ail 120011) . in which a red- 
dened, young starburst population can look similar to an older, 
unextincted stellar population, is also probably contributing 
to the ambiguity in the stellar ages. However, we expect that 
much of the difficulty in constraining stellar population ages is 
likely to arise from the application of relatively simple stellar 
models to what is likely an intrinsically complex mix of stel- 
lar populations of differing metallicity, reddening, and age. 
Accordingly, we could make use of stellar population models 
which employ more complex star formation histories, such 
as short bu rsts on top of an exponential star formation his - 
tory (e.g., iBorch et alj|2006t iKavirai et alj|2007t |Pve1l2008l) : 
however, we feel that doing so is not justified because the ob- 
servational data are insufficient to constrain even simple star 
formation histories and cannot provide a unique solution. We 
also stress that the uncertainties here, for SMGs with known 
redshifts, are far less than if the redshift also has to be derived 
from the photometric fit. 

The star formation history, and by extension, age, of a stel- 
lar population are critical elements in determining the ap- 
propriate M/L ratio needed to accurately convert a single 
galaxy's stellar luminosity to stellar mass. However, if one is 
willing to forego accurate masses for individual galaxies, one 
may take advantage of statistics to obtain reasonably precise 
average and typical stellar masses from a large sample of indi- 
vidual galaxy masses which are much more uncertain. Fortu- 
nately, the light-to-mass (L/M) ratios of stellar populations of 
different star formation histories at ages above 50 Myr, shown 
in Figure|2]for a sample of star formation histories, have a nar- 
rower range of values than the ages themselves and are thus 
less uncertain than the ages: for example, in the M05 stellar 
population models, the L/M ratio of an instantaneous burst 
(IB hereafter) population with age 100 Myr differs by only 
20% from that of a constant star formation (CSF hereafter) 
population with age 2 Gyr. At any age in the range 50 Myr 
to 3 Gyr (approximately the age of the universe at z ~ 2), the 
L/M of the star formation history with typically the lowest 
L/M ratio (IB) differs by only a factor of 2-3 at most from 
that of the star formation history with the highest typical L/M 
(CSF). So, with no constraint on the star formation history, we 
can still constrain the stellar masses of individual galaxies to 
within a factor of 2-3, neglecting other systematic uncertain- 
ties. If we assume the star formation histories of individual 
SMGs are uniformly distributed over the range of possible 
histories, which is reasonable because our attempts to fit the 
SEDs of SMGs with a variety of star formation histories re- 
vealed no single preferred history over the entire sample, the 
average star formation history for the population and the aver- 
age L/M ratio should exhibit less dispersion than those asso- 

13 HHik et al. (2003) show, using their sample of z ~ 2 massive red galax- 
ies with spectroscopic redshifts, that the best-fit ages obtained through stellar 
population synthesis modeling of SEDs constructed using broadband JHK 
photometry have a mean absolute deviation of a factor of ~ 3 relative to the 
ages derived from SEDs constructed using low-resolution near-IR spectra. 
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Figure 2. Lfj /M* as a function of time since onset of star formation for solar- 
metallicity population synthesis models of BC03 and M05. Solid lines rep- 
resent models of BC03, while dotted lines represent the models of Maraston 
1 2005). The increase in L///M* near log age = 8.5 in the Maraston ( 2005) 
burst model marks the onset of the TP-AGB phase. Because at any age in 
between 50 Myr and 3 Gyr the Lh /M* of the star formation history with typ- 
ically the lowest L#/M* ratio (IB) differs from that of the star formation 
history with the highest typical Lh/M„ (CSF) by at most a factor of 2-3, we 
can constrain the stellar masses of individual SMGs to within a factor of 2-3 
without knowing the star formation history. 

ciated with individual galaxies. Consequently, the uncertainty 
in the sample-averaged stellar mass derived using the average 
L/M ratio will be much less than the factor of 2-3 uncertainty 
in the individual masses. 

In the absence of any indication of which star formation his- 
tory is most appropriate, and with the knowledge that (a) the 
SED fits of the different star formation histories to our SMG 
sample are correlated; (b) the best-fit ages are likely to not be 
representative of the true stellar population since the light we 
receive is strongly modified by reddening; and (c) the L/M 
ratios implied by the different star formation histories will be 
similar to within a factor of a few, we elect to compute stel- 
lar masses for our sample which are averages over different 
possible star formation histories, instead of basing our com- 
putation on a particular star fo rmation history. We adopt the 
approach of Bo rvs et all d2005l) . in which we assume a single 
L/M ratio for the entire SMG sample for each of the IB and 
CSF histories, and average the resulting masses. We choose 
these two because they are expected to represent the lower and 
upper extremes of the stellar properties predicted by the dif- 
ferent star formation histories; in averaging their results, we 
hope to reduce the systematic effects of star formation his- 
tory so that our stellar mass estimates are correct to within a 
factor of 2. The mean and median of the stellar masses cal- 
culated in this way of the full SMG sample will thus be the 
best-constrained quantity we obtain. 

3. 1.2. Near-IR Excess in SMGs and Its Removal 

In Paper I, we found that numerous galaxies in our SMG 
sample display evidence of non-stellar emission in their rest- 
frame J — H and H-K colors (see Figure 6 in Paper I), by 
comparing to samples of "normal", non-active local galax- 
ies and local ULIRGs. Upon close examination of the ob- 
served SEDs of our sample, we find that approximately 20% 
of the sample display a clear red upturn in the rest-frame 
near-IR portion of their SED. The red upturn suggests the 
presence of hot dust and possibly an obscured AGN, since 
stellar SEDs typically fall past the 1 .6 /im stellar peak. An 
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additional 50% of the SMGs in the sample display smaller 
near-IR excesses redward of the 1 .6 /jm stellar peak similar to 
those found in 24 /mi-bright luminous IR galaxies and AT-band 
sele cted galaxies at 0.5 < z < 2.0 by Mag nelli et all (12008b 
and Mentu ch et alj Q2009). These authors suggest that the ex- 
cess arises from very small, hot dust grains in star formation 
regions, compar able to the near-IR emission detected from 
Galact ic cirrus (|Flagev et all 120061) and reflection nebulae 
(e.g.. lSellgren et al]|1983l) as well as local star-forming galax- 
ies (e.g.. iLu et al] 120031) . which can be described by gray- 
body emission at temperatures ~ 1000 K. 56% of our s am- 
ple also meet the criterion of lAlonso-Herrero et alJ (120061) for 
IRAC power-law galaxies (PLGs; power law slope a > 0.5), 
which those authors consider to be obscured AGNs. Yet, it is 
not completely clear that th e PLG classification can generally 
be attributed to AGN since iDonley etalJ (120081) find that star 
forming ULIRG templates at z > 2.6 meet the PLG criteria as 
well (indeed, three of the four SMGs that show no evidence 
of SED upturn but have a > 0.5 are at z > 2.8). Regardless 
of whether or not the origin of the observed excess continuum 
emission is an obscured AGN or star formation regions, its 
contribution to the near-IR emission of ~ 70% of the SMGs 
in our sample must be removed in order to derive the most 
representative stellar luminosity and mass. 

We remove the non-stellar continuum contributions to the 
near-IR emission of our SMG sample by assuming that the 
total observed SED can be represented as a sum of a stellar 
population model and a power law (PL) of the form f v oc v~ a 
representing the excess, where /„ is the flux density at fre- 
quency v. Note that, if the excess is in fact due to graybody 
dust emission in star forming regions rather than dust heated 
by an obscured AGN, the PL model will still be approxi- 
mately correct, since the modified Planck function describing 
a graybody can be approximated as a PL over a small range 
of wavelength (e.g., a graybody of temperature 700-1000K 
and emissivity t v oc v , appropriate for small dust grains in 
star-forming regions, can be approximated by a PL with in- 
dex ranging from -1.8 to -4 in the wavelength range 2-3 /im, 
depending on the temperature). However, if the excess is due 
to a 3.3 /im polycyclic aromatic hydrocarbon (PAH) emission 
line, the PL model will not improve the fits. We scale the 
PL to the IRAC 8.0 /im data point, and denote this scaled PL 
as the maximum possible PL contribution. We then subtract 
from all the observed photometry increasing fractions of the 
flux in this scaled PL, in steps of 0. 1 x (maximum PL contri- 
bution) and fit the BC03 and M05 stellar population models to 
the PL-subtracted data points. The PL fraction which results 
in the lowest value of xt gives us the luminosity of stars only, 
which we can then use to find the stellar mass. We estimate 
the error in the PL fraction by determining the PL fractions 
which result in a change of A\l = 1 above and below the 
best-fit PL fraction. 

The choice of the best PL index to use in our PL re- 
moval is not obvious, because differing obscuration lev- 
els can cause the rest-frame near-IR PL slope of hot dust 
emission to vary between gal axies even when the AGN 
powers the dust emission (e.g., Simpson & Raw lingsll2000l) . 
Observationally-derived near-IR PL slopes for local unob- 
scure d, type-1 AGN have a typ i cal value of a ~ 1 — 1 .5 
(e.g.. iMalkan & Filippenkol 119831: lEdelson & Malkanl 119861: 
Neugebaue r et alll987l) . whereas the range for local Se yfert 2 
nuclei is a ~ 2-4 (e.g.. lAlonso-Herrero et al.1 12003). Ob- 
scured starbursts will have similar near- and mid-IR PL slopes 



to obscured (type 2) AGN: as mentioned in the previous para- 
graph, graybody dust emission from star forming regions can 
be approximated over A ~ 2 - 3 /im as a PL with a ~ 2 - 4. 
We initially chose to use a PL index of a = 2, based on the 
observed average mid-IR continuum sl ope for SMGs with 
IRS spectra, «mir = 1 .97 ±0.22 (Men endez-Delmestre et alJ 
2009). However, for some galaxies in our sample, the v 2 PL 
is too shallow and results in a very poor fit; thus, we have 
allowed the PL in our continuum removal to have the form 
v~ 2 or v~^, whichever best fits the data in a \ 2 sense for each 
SMG. By allowing 2 different possible values of the PL index, 
we introduce a additional uncertainty into our stellar masses; 
howe ver, the uncertainty in the average mass is small (see 
33.3b . The results of our SED fits suggest that a = 3 typi- 
cally produces a better fit to the observed data overall: for the 
82% (72%) of the SMG sample which have a non-zero best-fit 
PL contribution at observed-frame 8.0/im using the M05 IB 
(CSF) models, 70% (79%) have best-fit PL indices of a = 3 
while 30% (21%) have best-fit PL indices of a = 2. When 
the BC03 stellar models are used in the SED fits, the relative 
fractions change somewhat: for the 70% (79%) of the sample 
with a non-zero best-fit PL fraction for the IB (CSF) star for- 
mation history, 74% (8 1 %) have best-fit PL indices of a = 3 
while 26% (19%) have best-fit PL indices of a = 2. 

We note that when we extrapolate the observed-frame 
24 /im flux from the 8.0 /im flux based on the best-fit PL frac- 
tion, the extrapolated value exceeds the observed 24 /im flux 
in ~ 50% of cases when a = 3 (the median over-prediction for 
that 50% of the sample is a factor of 2.3) and for ~ 20% of the 
sample when a = 2 (by a median factor of 1.5). The ~ 50% 
and ~ 80% which are not over-predicted by a = 3 and a = 2 
PLs, respectively, are under-predicted. The over-prediction of 
the rest-frame mid-IR flux based on the rest-frame near-IR PL 
implies that the near-IR and the mid-IR require different PL 
slopes. This phenomenon i s observed in local Seyfert ga laxies 
(e.g., lEdelson et all 119871: lAlonso-Herrero et al] 120031) . and 
can occur when an unobscured nucleus' emission is modi- 
fied by a foreground reddening screen (e.g., its host galaxy). 
Such cur vature in the m i d-IR SEP of high-z QSOs is also 
found by Galla gher et alJ d2007l) . who find that the rest-frame 
8.0 /im luminosities of their sample are over-predicted by fac- 
tors ranging from 1.25-2.5 upon extrapolating the best-fit 
observed-frame 1.8-8.0/im PL to rest-frame 8.0 /im. Thus, 
our choice of near-IR PL indices for fitting are not necessarily 
invalidated by over-predicting the observed 24 /im flux. 

An alternate approach to fitting the near-IR excess emis- 
sion in the observed SED for the purpose of estimating stellar 
masses is to fit a combination of stellar population synthesis 
models and AGN (type 1 and/or 2) t emplates from the lit- 
erature (see, e.g.. iMerloni et al]|2010l) . However, we prefer 
the stellar population plus PL approach over the use of tem- 
plates for several important reasons. First, because we have 
no spatial information for the near-IR excess observed in our 
SMGs, and because SMGs are strongly star-forming galaxies, 
we are hesitant to assume a priori that the excess arises from 
an AGN when dust heated by star formation could also ex- 
plain the observations. The vast majority of our SMG sample 
are adequately fit in the rest-frame UV and optical by stel- 
lar models only, so it is not unreasonable to consider the case 
in which stellar radiation powers the near-IR excess. Setting 
aside our wish to not make assumptions regarding the origin 
of the near-IR excess, our observation that the SMG SEDs 
can generally be described by stellar emission argues against 
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a need for type-1 AGN templates; a case could still be made, 
though, for the use of type-2 AGN templates. On the other 
hand, our stellar population model plus PL representation is a 
good reproduction of a type-2 AGN SED, since a type-2 AGN 
SED is dominated by stellar light in the rest-frame optical and 
PL continuum in the near-IR. Moreover, by using our stellar 
model plus PL approach it is actually more straightforward to 
separate the contributions of stellar emission and featureless 
continuum to an SED than would be possible using a stellar 
model and a type-2 AGN template, because the type-2 AGN 
template may contain some contribution from the AGN host 
galaxy in the near-IR. Finally, and most importantly, we prefer 
to use our PL removal model rather than AGN templates be- 
cause the PL model permits some flexibility in the near-IR PL 
index. We discussed above how our efforts to determine the 
most appropriate near-IR PL index for the removal of the con- 
tinuum excess revealed that we could not fit all of the SMGs 
in our sample with a single PL index; PL index variation in the 
near-IR is also observed in local AGN. Yet, by using an AGN 
template to fit the observed SEDs, we would have only one 
option for the near-IR PL index, which may or may not be ap- 
propriate for our sample. For example, the average optically- 
bright (type-1) QSO SED constructed f rom Sloan Digital Sky 
Survey QSOs by Richard s et al.l ( 120061) . frequently used as an 
AGN template, has a PL index a = 1 .8 in the wavelength range 
1-3 fim. This PL index will somewhat under-subtract the 
near-IR continuum component in those SMGs which are best 
fit by a = 2 and significantly under-subtract the continuum 
for the SMGs which are best-fit by a = 3 (the majority of the 
sample ). The composite type-2 AGN SED of Pollet ta et al.l 
(2007), another commonly-used AGN template, has a near- 
IR PL index a = 1.3 between 1 and 3 fim; this template will 
thus significantly under-subtract the continuum component in 
nearly our entire SMG sample. 

We illustrate the effects of our PL removal procedure on a 
typical SMG in our sample showing a near-IR upturn in Fig- 
ure [3] We first show the observed optical-IRAC photometry 
and the best-fit CSF stellar population model from M05 to the 
observed data; the xt of the fit indicates that the fit is poor. 
We then show the best-fit combination of v~ 2 PL and CSF 
stellar population model to the same observed photometry as 
well as the flux contributed by the PL component. Finally, 
we show only the stellar component of the fit from the middle 
panel which allows us to determine the stellar luminosity of 
the SMG: the PL-subtracted photometry is plotted along with 
its corresponding best-fit CSF model. Note the significant im- 
provement in the xt statistic of the stellar model fit when the 
PL component is subtracted as well as the more appropriate 
age of the best-fit model. For ~ 60% of the SMGs the best xl 
value of the fit improves when we subtract a PL contribution 
from the observed SED; indeed, for ~ 20% of the sample, the 
best xl °f me fit decreases by a factor of 2 or more when we 
subtract the PL contribution. In Figure|4]we show the results 
of the best CSF stellar population plus PL fit to the observed 
photometry for the SMGs in our sample requiring PL contri- 
butions at observed-frame 8.0 /im greater than 70%, for which 
the fits improved the most with the inclusion of the PL compo- 
nent. In Figure|6]we show the best CSF stellar model plus PL 
fits for the remainder of the SMGs, for which the necessity of 
a PL SED component can be more ambiguous. We note that, 
for both Figures [4] and [6j the IB fits are visually similar to 
the CSF fits, although the best-fit PL contributions differ be- 
tween the IB and CSF star formation histories for ~ 60% of 
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Figure 3. Example of stellar population model plus power-law fraction 
SED fitting used for our sample of SMGs. Panel (a): The observed 
UBVRIz'JK+lRAC photometry (round points) for a typical SMG with a near- 
IR excess is shown with the best-fit CSF stellar population model of M05 to 
that data. The near-IR excess is most clearly seen at \ rest > 2fim. We note 
that the requirement that the age of the best-fit model be physical was relaxed 
in this particular case to allow the fit to converge. Panel (b): The observed 
photometry for the same SMG is shown along with the best-fit combination 
of i/~ 2 power law and CSF stellar population model from M05. Panel (c): The 
power-law subtracted photometry (black stars), representing the total stellar 
contribution to the observed SED is shown with the best-fit CSF stellar popu- 
lation model from M05 to those data points only. Note from panel (a) that the 
unphysical age of the best-fit stellar population model (it is greater than the 
age of the universe at z = 2.4) and the statistically poor fit are indications that 
the stellar population model alone is insufficient to account for the near-IR 
emission from this SMG. 

our sample for both the BC03 and M05 models. In contrast to 
the similar visual appearance of the IB and CSF fits, the best 
xt values for the IB star formation history are 10% greater 
(BC03 models) and 16% greater (M05 models), on average, 
than the best xt values for the CSF star formation history. 

3.2. Stellar Mass Calculation Details 

We calculate a stellar mass (M + ) for each SMG from the 
de-reddened, best-fit-power-law-subtracted absolute H mag- 
nitude (Mh). /f-band light has been used to determine the 
stellar mass in past studies of the s tellar mass of SMGs (e.g., 
iBorys et alJ200 5: Michalowski et al. 2010) because of its low 
sensitivity to previous star formation history and dust obscu- 
ration. We instead choose to use the rest-frame H band to es- 
timate M+ for our sample of SMGs to benefit from the low ex- 
tinction of the near-IR bands while minimizing the uncertain 
influence of thermally pulsating asymptotic giant branch (TP- 
AGB) stars on the model /f-band light-to-mass (L/M) ratio 
in the M05 stellar population models. The change in //-band 
L/M ratio induced by the appearance of a dominant solar- 
metallicity TP-AGB population at ages of ~ 0.5- 1.5Gyr is 
a factor of ~ 2, smaller than the change in the Zf-band ra- 
tio. Calculating the mass from Mh also has the advantage that 
the contribution of the power-law emission is naturally lower 
in the //-band than in /if -band, so the uncertainty introduced 
due to our subtraction of the power-law component is smaller. 
A final reason for the use of //-band in deriving the stellar 
masses is that for the highest-redshift (z > 3) SMGs in our 
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Figure 4. Rest-frame UV through near-IR SEDs of spectroscopically-identified SMGs in our sample unambiguously requiring significant near-IR power law 
components to fit the observed SED. Round points are observed photometric data, while non-detections are indicated by downward arrows originating from the 
measured upper limit. The best-fit CSF model of M05 is over-plotted as the black line. The best-fit combined stellar population model plus best-fit power law is 
represented by the dotted (blue/gray) line. The UV excess displayed by several of these SMGs cannot be reproduced by stellar or near-IR power-law components 
or a combination of the two. 
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Figure 5. Figure f4]continued. 
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Figure 6. Rest-frame UV through near-IR SEDs of remaining SMGs in our sample. Symbols and lines are as in Figure \4\ The large variation in non-stellar 
near-IR continuum contribution within the sample is apparent when comparing to Figure|4] ranging from pure stellar emission to nearly pure power-law emission. 



sample, the longest-wavelength IRAC observations (8.0 /im) 
do not quite sample rest-frame K, so calculations of the abso- 
lute K magnitude (Mk) for these sources are not constrained 
by th e avail able data without extrapolation. 

In {33.1. II we explained o ur rationale for adop ting the stellar 
mass calculation method of iBorvs et alj d2005l) . which uses a 
single L/M ratio for the entire sample of SMGs for a given 
star formation history. We obtain this Lh/M* ratio for each 
star formation history (IB and CSF) by finding the average of 
the best-fit model ages over the SMG sample for the fits of that 
star formation history and taking the L/M ratio for that star 
formation history corresponding to that average age. Thus, for 
the BC03 models we use L H /M± = 5.8L Q M^ for the CSF his- 
tory and Lh/M* = 7.9LqMq for the IB history, correspond- 
ing to average sample ages of 1.5 Gyr and 140 Myr, respec- 
tively. From the M05 stellar models we derive an average 



age and L jM ratio of 1 .6 Gyr and L H /M+ = 7 . 9 L© M Q ' for the 
CSF history, and 120 Myr and L H /M+ = 8.5 LqM" 1 for the IB 
history. 

We calculate the stellar mass for each SMG for each of the 
CSF and IB star formation histories using the Lh /M+ ratio for 
that history and M# from the normalization of the SED fit of 
that star formation history, then average the masses obtained 
from the IB and CSF histories to obtain our final stellar mass 
estimate for each SMG. We estimate the uncertainty in the 
stellar masses by dividing in half the difference in the masses 
resulting from the CSF and IB star formation histories, since 
these values tend to represent the maximum and minimum 
values of the Lh/M* ratio among the various simple star for- 
mation histories we examined (see Figure|2]l. We thus attempt 
to represent the systematic effects of star formation history 
and Lh/M* ratio in our uncertainty estimates, which are much 
more significant than the random errors in the model fitting 
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Figure 7. Figure [6]continued. 
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due to the photometric uncertainties (see, e.g., iMuzzin et alJ 
2009). However, we caution that myriad systematic errors are 
still possible for our mass estimates based on the assumptions 
we have made, as we discuss in the next section. 

3.3. Sources of Systematic Uncertainty in Median M+ 

Although the median stellar mass of our SMG sample is the 
best-constrained quantity we derive in our study, the median 
stellar mass still has associated random and systematic uncer- 
tainties which can be significant and must be mentioned. We 
have estimated the random error in our median stellar mass 
by randomly perturbing the observed photometry of our sam- 
ple within its associated error bars and carrying out our SED 
fitting and mass estimation procedure on the perturbed data 
for the galaxies. By analyzing the statistics over many tri- 
als (where each trial produces a median stellar mass for the 
sample), we find that the standard deviation of the median 



stellar mass, which we call the random error in the median 
of our sample, is 11%. Because the random error associated 
with our SED fitting is only ~ 10%, the overall uncertainty 
in our median stellar mass is dominated by systematic un- 
certainty. On the whole, we expect the dominant systematic 
uncertainty arises from our use of stellar population synthe- 
sis models to obtain a L/M ratio; as this uncertainty is dif- 
ficult to quantify and affects nearly all studies of the stel- 
lar mass in high-redshift galaxies similarly, we refrain from 
discussing it at length. Rather, we restrict ourselves to a 
brief summary of the more concrete systematic effects and 
their impact on our median stellar mass, noting that many 
of the sources of uncertainty have been explored in detail 
elsewhere (e.g..|Papovich et alJ200illvan der Wel et al.f2 006: 
IConrov et alj|2009tlMuzzin et al.ll2009l) . We further note that, 
in principle, systematic errors in the stellar mass calculations 
could occur within our sample between the different sky fields 
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because some of the fields have fewer optical photometric 
data points to constrain the SED models from which derive 
the ages and M/L ratios. However, in practice, since our aim 
is to constrain the average stellar mass over the SMG sam- 
ple and not individual galaxy masses, this particular source of 
systematic error is not significant: tests varying the number of 
optical photometric data points used in fitting the GOODS-N 
subsample of SMGs (which has the most photometric data of 
all the fields in our sample) reveal variation of no more than 
0. 1 dex in the mean and median stellar mass. 

Relatively minor contributions (~ 10-20% individually) 
to the systematic uncertainty in our median stellar mass arise 
from our metallicity and extinction law assumptions and our 
removal of the near-IR continuum exce ss in fitting the SEDs 
of our sample of SMGs . As shown by iMuzzin etafl (120091) 
and lConrov et al.l (120091) . the assumption of a single metallic- 
ity in our fitting introduces uncertainties in the stellar masses 
of 10-20%. We have estimated the uncertainty in the me dian 
stellar mass resulting from use of the lCalzetti et al.l ( 120001) ex- 
tinction law for starbursts to be also 10-20% by fitting stellar 
population models to the SEDs of our sample assuming the 
extinction laws for the Mil ky Way (Fitzpatrick 1986) a nd the 
Smal l Magellanic Cloud dPrevot et al.l 11984 iBouchet et all 
1985). We note that our simplistic treatment of the extinc- 
tion within SMGs as a uniform dust screen also affects our 
stellar mass estimates in a way which is difficult to quantify; 



yet, as there are essentially no constraints on the distribution 
of dust and star formation within SMGs, it is difficult to jus- 
tify the use of more compl ex extinction models, s uch as age- 
dependent extinction (e.g., Poggia nti & Wul l2000). The sys- 
tematic uncertainty in our average stellar mass estimate asso- 
ciated with our near-IR continuum removal procedure arises 
from the index of the PL subtracted from the observed SED 
to obtain the stellar H-band luminosity. We have estimated 
this uncertainty as 10-15% by comparing the stellar masses 
of the individual galaxies in our SMG sample resulting from 
the assumption of PL indices a = 2 and a = 3. 

The systematic uncertainty in our typical stellar mass esti- 
mate from our choice of IMF is considerably larger than those 
associated with metallicity, extinction, and choice of PL in- 
dex. The choice of a particular function effectively introduces 
a scaling factor into the mass estimates: for example, if we in- 
stead chose to use the Salpeter IMF wit h the BC03 and M05 
stellar populat ion models rather than the Chabrier d2003l) and 
Rroupa ( 2001 ) IMFs, the median stellar mass would be higher 
by a factor of 1.6 (Rroupa) and 1.8 (Chabrier). If we assume 
the systematic uncertainties arising from the choice of IMF, 
metallicity, extinction law, and PL subtraction are all indepen- 
dent (which may not be correct) and add them in quadrature, 
we obtain a conservative estimate of the systematic uncer- 
tainty in our typical stellar mass of a factor of 1.6-1.8, which 
is clearly dominated by the IMF systematic. Including the un- 
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certainty due to the stellar population models will obviously 
increase the overall systematic error estimate. 

We caution that a final source of systematic uncertainty 
which is not easy to quantify with currently existing obser- 
vational data may significantly affect our stellar mass esti- 
mates. Should a large fraction of the stellar content of SMGs 
lie within or behind the heavily obscured star formation re- 
gions, we are unable to trace its mass since we do not receive 
its light. In principle, dynamical mass estimates can help as- 
sess the amount of mass missing from our estimates, but the 
currently available dynamical masses for SMGs come from 
CO or Ha observations which may not trace the bulk of the 
stellar mass. In the future, when precise dynamical masses for 
our sample are obtained from rest-frame near-IR (i.e., stel- 
lar) light, we will be able to fully address this issue. Until 
then, our stellar masses cannot include this unknowable, but 
potentially significant, fraction of the stellar mass in SMGs. 
How much mass is missing from our estimate of the median 
will depend on how the dust is distributed throughout SMGs, 
which we discuss further in the next section. 

4. RESULTS 

In Table |2j we provide for our sample of 65 
spectroscopically-identified SMGs the de-red dene d val- 
ues of Mh and the calculated as described in §13.21 as well 
as averaged V-band extinctions and the 8.0 /im PL fraction 
resulting from our stellar population synthesis model plus PL 
fitting for the M05 stellar models. The results for the BC03 
models are largely similar, though with a scaling factor in 
the stellar masses, so we have chosen not to tabulate those 
results here. As a brief summary of the differences between 
the BC03 and M05 results, we find that the absolute //-band 
magnitudes resulting from the use of the BC03 models are 
0.03 mag fainter, on average, than those resulting from the 
M05 models; the V-band extinctions are equal, on average, 
between the BC03 and M05 models for the IB star formation 
history but the extinctions for the CSF star formation history 
are 0.4 mag higher systematically in the BC03 results; 
finally, the stellar masses resulting from the BC03 models 
are systematically 25% higher than those obtained with the 
M05 models, although we note that a Chabrier IMF has been 
used for the BC03 results while a Rroupa IMF has been 
used for the M05 results. We discuss the results contained 
in Table [2] some of their implications, and some caveats in 
the sections below, providing the quantities derived from 
the BC03 models only when significantly different. For our 
discussion we separate the z < 1 .5 SMGs from the higher-z 
galaxies because our focus in this paper is on the high-redshift 
SMGs which dominate our sample, cover a smaller epoch 
of cosmic time, and are more important in a cosmological 
sense in that they contribute more significantly to the cosmic 
star formation rate density (z ~ 2 SMGs with 5s5o > 4mJy 
contribute ~ 25% of the overall star formation rate density at 
that redshift, while similarly bright z < 1 SMGs contribute 
~ 10%; C05; Wardlow et al. 2011). Moreover, the low-z 
SMGs are known to differ somewhat from the high-z SMGs 
in that they are typically less far-IR-luminous (C05), and thus 
may have different characteristics from the high-z sub-sample 
overall. 

4.1. Contribution of Non-Stellar Near-IR Emission in SMGs 

Through fitting the combination of stellar population mod- 
els and a PL to the observed SED of our sample of SMGs, we 
have placed constraints on the fraction of the observed-frame 
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Figure 9. Cumulative distribution of the fraction of rest-frame //-band lumi- 
nosity in the power-law component for z > 1.5 SMGs, based on the best-fit 
power law fraction at observed-frame 8.0 (im, using the M05 stellar popula- 
tion models. Approximately half the SMG sample has less than 10% con- 
tribution to the rest-frame H luminosity from the power-law component, and 
11% of the sample has more than 50% of the //-band emission originating 
from the power-law component component. 



8.0 /im emission of each SMG which is contributed by a PL 
component, which we tabulate in Table|2] To make our results 
easier to interpret, we have converted this best-fit PL fraction 
at observed-frame 8.0 /im into the fraction of the total rest- 
frame //-band luminosity (L#; approximately the location of 
the peak of the stellar emission) of each SMG in our sample 
contributed by its PL component. In Figure|9]we show the cu- 
mulative distribution of the calculated fraction of rest-frame 
H luminosity in the PL component for the z > 1 5 portion of 
our sample for both the IB and CSF star formation histories 
using the M05 stellar population models. The distributions of 
inferred PL fraction are quite similar for the different star for- 
mation histories; for each, approximately half of our z > 1.5 
sample have essentially negligible contributions to their total 
Lh from a PL component. Of our entire spectroscopic SMG 
sample, 69% (79%) of the galaxies have non-zero fractions 
of Lh for the CSF (IB) star formation history; yet for only 
1 1 % (for both star formation histories) is the fraction of Lh in 
the PL component greater than 0.5. For reference, and warn- 
ing, we list in Table [3] those SMGs for which the fraction of 
Lh in the PL component is greater than 0.5. Clearly, the ma- 
jority of our SMG sample are stellar-dominated in rest-frame 
H, although for much of the sample the stellar mass will be 
somewhat overestimated (by a median value of 10%, and a 
mean value of 25%) if we assume that all of the near-IR light 
is stellar. While the amount by which the stellar masses are 
overe stim ated lies well within the errors of our present study 
(see {33.31 l. future studies which aim to determine more 
precisely for individual SMGs will need to remove the non- 
stellar contribution to the near-IR light. 

Interestingly, one of the near-IR PL-dominated SMGs 
listed in Table [3] also displays a rest-frame UV excess. 
At least three additional galaxies in our sample also ap- 
pear to have a UV excess, which are noted in Table |2 
we note that because the detection of a UV excess re- 
quires a reported observation in the observed-frame U band, 
which most of sample lack, we are prevented from accu- 
rately determining the frequency with which UV excesses 
occur in our sample. Each of the four UV-excess SMGs 
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have significant PL contributions in the near-IR, even though 
three of the four galaxies do not qualify as PL-dominated. 
For two of these SMGs (SMM J123632.61+620800.1 and 
S MM J123635. 59+621424.1) AGN spectral i dentifications in 
the ob served-frame optical (C05) and X-ray (Ale xander et alj 
120051) suggest that a type-1 AGN could be the origin of 
the UV excess; high-resolution imaging of these galaxies in 
the observed-frame //-band showing a compact morphology 
would lend additional support to such a hypothesis. However, 
the other two UV-excess SMGs (SMM J123553. 26+621337.7 
and SMM J 123555. 14+62090 1.7) have no AGN signatures in 
their observed-frame optical spectra, so attributing the UV ex- 
cess to a type-1 AGN is problematic. For these SMGs, a pos- 
sible explanation for the UV excess is that a simple uniform 
dust screen is simply an inappropria te model of the stellar 
extinction across the entire galaxy, as Pogg ianti et al.l (2001) 
show for luminous IR galaxies displaying e(a) spectra. For 
example, in these particular SMGs some portion of the total 
star formation may be relatively unobscured and contribute 
some UV flux to the SED, while the dominant star formation 
activity is heavily obscured. Another possibility is that the 
UV excess is a result of scattering of obscured AGN light, 
similar to the origin of the far-UV flux in the local Seyfert 2 
galaxy NGC 1068 (e.g. JCode et al.ll993tlGrimes et al.ll999h . 
High-resolution U -band imaging of these SMGs will also be 
important to determine the origin of the UV excess. 

4.2. Near-IR Stellar Luminosity of SMGs 

As a result of our separation of the PL component from 
the stellar component in the observed SEDs of our sample of 
SMGs, we have obtained more representative estimates of the 
rest-frame near-IR stellar luminosity for SMGs. The abso- 
lute H magnitudes resulting from our SED fitting procedure 
are much less dependent on star formation history than the 
stellar masses, as the absolute magnitudes are calculated from 
the photometric data point nearest the rest-frame H band by 
applying a (small) /f-correction interpolated from the best-fit 
stellar population model. Thus, the absolute magnitudes of in- 
dividual SMGs listed in Table [2] are much better constrained 
than the stellar masses listed for individual galaxies. Since we 
fitted both IB and CSF star formation histories to the observed 
SMG SEDs, we obtain multiple estimates of Mh', to obtain a 
single Mh for each SMG for our present discussion and for 
Table [2] we average the Mh values resulting from the fits of 
the IB and CSF star formation histories. We estimate the un- 
certainty as half of the difference between the IB and CSF 
values; by doing so we take into account the uncertainty in 
the .^-correction determined from the SED models of differ- 
ent star formation history, although we anticipate it is small. 

The median of the stellar absolute H magnitude determined 
in this way for the z > 1.5 SMGs in our sample, including 
the upper limits for galaxies which had fewer than 3 detec- 
tions in their SEDs, is (Mh) = -25.7 ±0.2 with standard devi- 
ation om h = 1.1 for the M05 models; we note that the median 
Mh found using the BC03 models is only 0.02 mag fainter. 
When de-reddened by the median value of Ay, the median ab- 
solute //-band magnitude for becomes (Mh) = -26.0. For the 
z < 1 .5 SMGs, we find a median, uncorrected for reddening, 
of (Mh) = -25.3 ±0.7 with standard deviation om h = 1.7 for 
the M05 models, and when de-reddened by the median Ay, the 
median becomes (Mh) = -25.6. There is thus no statistically 
significant difference in the typical Mh between the low- and 
high-z sub-samples. Both the low- and high-z sub-samples of 
SMGs have typical stellar luminosities which are greater than 



local L* galaxies (M H = -23.7; ICole et al.ll2001l) , which sug- 
gests that the stellar populations of SMGs are younger, that 
typical SMGs are more massive than local L* galaxies, or 
some combination of the two. Given that we observe these 
SMGs at high redshift, it is highly likely that their stellar pop- 
ulations are younger than those of old local L* galaxies, but 
we should compare the stellar masses of SMGs to those of L* 
galaxies directly to examine the effect of mass. 

4.3. Stellar Mass Results for SMGs 

As mentioned in $3.21 we have chosen a stellar mass esti- 
mation procedure which constrains the average over the entire 
sample at the expense of accuracy in the mass of the individ- 
ual galaxies. As a result, although we list stellar masses for 
individual SMGs in Table |2]for readers who w ish t o use them 
in spite of their uncertainty (factor of ~ 2; see 93.2t . we do not 
wish to focus on results in an individual sense. We find it more 
appropriate to discuss the average (median/mean) mass of the 
sample and its range, which are not as sensitive to the star 
formation history uncertainty due to our large sample size. 
In Figure [10] we show the distribution of for the z > 1 .5 
SMGs in our sample for both the M05 and the BC03 mod- 
els to visually display the range in stellar mass of our sample. 
We caution that we have included in the histograms the stellar 
masses even for SMGs whose rest-frame optical and/or near- 
IR flux appeared to be dominated by a PL component prior to 
its removal, and remind the reader th at the M05 ste llar mass 
estimates use a different IMF (that of Kroupa 2001) from the 
BC03 stellar mass estimates (that of Chabriel l2003l) . While 
the distributions for the M05 and BC03 models have similar 
widths, a ~ 0.4 dex, the distribution of BC03 stellar masses 
is shifted relative to the distribution of M05 masses such that, 
even when corrected to the same IMF, the M05 masses are 
systematically ~ 20% lower than the BC03 masses on aver- 
age. We attribute at least some of the difference in the masses 
estimated using the M05 and BC03 models to the effects on 
the Lh/M+ ratio of TP-AGB stars included in the models of 
M05. At the average age of SMGs in the CSF model for M05, 
1.6 Gyr, the L H /M+ ratio is ~ 40% higher than that for the 
BC03 models at the average SMG age, 1.5 Gyr. We adopt 
the masses from the M05 models for the remainder of our 
discussion since the inclusion of the TP-AGB phase in stel- 
lar evolution in those models makes them more likely to be 
representative of the actual L/M ratio. 

The median of the stellar mass distribution for the z > 1 .5 
SMG sample is (M*) = (7.2 ± 1.2) x 1O 1O M for the M05 
models [(M*) = (9.0 ± 1.7) x 10 10 M Q for the BC03 mod- 
els]. The inter-quartile range (25th - 75th percentile) is 
3.3 x 10 10 - 1.3 x 1O U M , suggesting that the z > 1.5 SMGs 
span a range of approximately an order of magnitude in stel- 
lar mass. The z < 1.5 SMGs in our sample seem slightly 
less massive on average than the high-z SMGs, having a me- 
dian stellar mass of (M*) = (4.4 ± 1.7) X 10 1() M Q for the M05 
models; however, the difference in median mass between the 
low-z and high-z SMGs is only ~ 2a, and so not statistically 
significant. We note that the new median we obtain for 
the z > 1 .5 SMGs lies between that o f model z ~ 2 SMGs in 
the semi-analytic model examined in Swin bank et al.1 (I2008L 
= 2. 1 x 10 1() M Q ) and the model SM Gs in the hydrodynam- 
ical simulations of lDave etaT] d2010l M t = 2.7 x 10 n M Q ); 
thus our new estimate does not favor one set of theoretical 
models over the other. We also note that the median stellar 
mass we obtain here forz > 1.5 SMGs does not qualitatively 
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Figure 10. Histograms of stellar mass (A/*) for C05 sample of SMGs using 
the stellar population synthesis models of M05 (top) and BC03 (bottom). The 
filled histogram indicates the M* distribution when no power-law component 
is subtracted as part of the SED fitting; in those distributions a slight shift to 
higher M* can be observed. The distributions span approximately an order of 
magnitude, and the power-law-subtracted distributions have medians (7.2 ± 
1.2) X 1O 1O M (M05 models) and (9.0 ± 1.7) X 10 10 M Q (BC03 models). 



change the conclusion of lAlexander et alJd200 8a) that the ma- 
jority of SMGs fall a factor of 3-5 below comparably massive 
local galaxies on the Mbh~M+ relation; thus, any lag in the 
growth of the central black holes relative to the host galaxy is 
modest in SMGs. 

The typical stellar mass of SMGs found in previous obser- 
vational studies span an order of magnit ude; our new result 
is intermediate between the estimate of ISmail et al.l d2004l 
M, ~3x 10 I0 M Q ), who lacked rest-frame near-IR data for 
their analysis and thus could only weakly constrain the stel- 
lar mass, and the results of previous studi es which used rest- 
frame near-IR data from Spitzer-IRAC (IBorvs et al.l 120051; 
iDye et al.ll200l IMichalowski et al.1 120101) . Although lower, 
our median M* is marginally consistent with that obtained 
by iBorvs et ail 



3.3 ± 1.2 x 10 n M Q ), whose 
SMG sample is contained fully within ours, falling just out- 
side their ± 2a uncertainty interval. A better comparison to 
IBorvs et afl is to compute the median for only the subsam- 
ple of SMGs which were also studied by those authors; when 
we do so, we obtain (M*) = (1.3 ±0.3) x 10 11 M P) , w hich falls 
within the ±2rj uncertainty interval of Borvs et al.. We note 
that since the stellar population models we have used fo r our 
st udy are not avai l able f or the IMF used by Bor vs et al.l that 
of lMiller & Scalol (119791) . we cannot make a precise compari- 
son between the two results; however, the smaller contribution 



of low mass stars in the iMiller & Scalol d 19791) IMF relative 
to the Kroupa IMF suggests that the discrepancy in the me- 
dia n between the two stud ies may be somewhat larger under 
the IMiller & Scalol (119791) IM F To compare our new resul t 
to those of IDye et al.l (120081) and IMichalowski et al l d2010l) . 
who use a Salpeter IMF, we have also carried out our SED 
fits with a Salpeter IMF and find that our median mass for 
Z > 1.5 SMGs increases to 1.3 x 10 n M Q . Despite the in- 
crease resulting from the IMF conversion, our typ ical stellar 
mass is still significantly smaller than that found bv lDve et al.l 
(120081. (M^) = 63t] j x 10"^^), falling outs ide their ±3a 
uncertainty interval. Michalo wski et al.l (120101) obtain a me- 
dian stellar mass of (M*) = 3.5 x 1O U M , a factor of 2.7 
larger than our median M±\ the sig nificance of the difference 
is unclear since Michalo wski et al.l provide no uncertainty es- 
timates for their individual or median stellar masses. 

While the fact that we have used different stellar popula- 
tion mode ls ("e.g.. IBorvs et al|2 005 use the stellar population 
models of Bruzual & Chariot 1993) likely plays a role in our 
systematically lower stellar mass estimates, the uncertainties 
associated with the different models are insufficient to fully 
account for the discrepancies. One important factor underly- 
ing our smaller stellar mass estimates is that we do not as- 
sume that the near-IR light from SMG s is en tirel y stellar in 
origin , unlike the studies of B orvs et al.l d2005l) and lDve et al.l 
(2008). Rather, we have subtracted the non-stellar near-IR ex- 
cess from the observed photometry and used M# to calculate 
the stellar mass instead of Mk since the contribution of the 
power-law emission is lower in the H-band than in K-bwd. 
If we instead estimate M t for our sample using rest-frame K- 
band absolute magnitudes and do not subtract the near-IR ex- 
cess, the median mass for our sample is (M*) = 1.3 x 1O U M 
under the M05 models and Kroupa IMF, a factor of ~ 2 
higher, reflecting the impact of the near-IR continuum excess 
on the stellar masses in an average sense. We thus expect 
that the previous studies in which the near-IR continuum ex- 
cess was not subtracted prior to stellar mass calculation have 
systematically overestimated the stellar masses of SMGs with 
near-IR excesses. 

In 34.21 above we noted that typical SMGs have brighter 
stellar Mh than local L* galaxies, which could be due to a 
higher stellar mass in SMGs or younger stellar population 
ages. It is thus notable that we find that the typical stellar mass 
of high-z SMGs is a factor of 2 below the typical stellar mass 
of local L* galaxies (m* - 1.4 x 10 n M ? ; ICoie et al.ll2001l) . 
Thus, the difference in stellar M H between SMGs and local 
L* galaxies can be attributed to the younger stellar population 
ages of SMGs. We note, though, that given the systematic 
error s in our mass estimates, which we discussed earlier in 
33.31 we consider the typical stellar mass of high-z SMGs to 
be roughly consistent with m*. We w ill re turn to the compar- 
ison of SMGs with local galaxies in 35.41 

4.4. Average Extinction of the Stellar Populations of SMGs 

Our stellar population model fitting, which assumes a sim- 
ple dust screen model, also provides estimates of the typi- 
cal V-band extinction for our sample of SMGs. Because the 
SMGs are unresolved in their observed-frame optical-IRAC 
photometry, we regard these extinctions as averages across 
the entire galaxy; the true obscuration probably varies with 
location. However, since the extinction estimates are corre- 
lated with stellar population age and star formation history, 
which we could not firmly constrain for individual SMGs, 
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Figure 11. Histograms of Ay for C05 sample of SMGs using the stellar popu- 
lation synthesis models of M05 for two different choices of star formation his- 
tory, assuming the Calzetti et al. (2000) extinction law for starburst galaxies. 
The total near-IR light we observe from SMGs is only modestly extincted, 
having median (Ay) = 1.2— 1.5, depending on the choice of star formation 
history. 



we once again consider the average and range over the en- 
tire sample to be the best-constrained quantities. We show the 
distribution of the fitted extinctions for both IB and CSF star 
formation histories for z > 1.5 SMGs in Figure [TT] which il- 
lustrate both the likely range of Ay over the sample and the 
effect of star formation history on the extinction estimates. 
The distributions for both star formation histories indicate 
that the galaxy-averaged optical and near-infrared light that 
we detect from high-z SMGs is typically moderately obscured 
(0.5 < Ay < 1.5), with medians of (Av.burst) = 1.5 ±0.1 mag 
and (Av.csf) = 1.2 ±0.2 mag and standard deviations of 0.8 
and 1.0 mag, respectively. The extinction in the z < 1.5 
sub-sample spans a similar range, with medians (Av.burst) = 
1.6 ±0.2 mag and (Av,csf) = 1.8 ±0.3 mag and standard de- 
viations of 0.8 and 0.7 mag. To give the reader a sense of the 
difference in Ay between star formation histories for individ- 
ual SMGs, the values we list in Table [2] are the average of 
the Ay's found for the CSF and IB star formation histories, 
and the errors are taken to be half the difference between the 
two. The typically modest visual extinctions we obtain are 
consistent with tho se estimated from the B aimer decrement 
dTakata et alj|2006l) . but contrast sharply with the extreme IR 
luminosities of SMGs, which imply powerful energy sources 
hidden behind copious dust, and with the typical visual ex- 
tinction derived from the optical depth of the 9.7 fj,m silicate 
absorption feature in the mid-IR spectra of SMGs (Ay ~ 6; 
Mene ndez-Delmestre et al.ll2.009t) . However, we may recon- 
cile the optical/near-IR Ay estimate with the mid-IR estimate 
by noting that the mid-IR light is less sensitive to dust ex- 
tinction so we receive emission from more obscured regions 
at mid-IR wavelengths than at optical wavelengths, thus in- 
creasing the average Ay we derive. 

The discrepancy between optical and mid-IR estimates of 
Ay highlights the difficulty in quantifying the stellar mass be- 
hind heavily obscured star formation complexes by means of 
an UV-near-IR SED fitting analysis, and again we caution that 
such mass may be missing from our estimates in this paper. 
How much the mass has been underestimated depends criti- 
cally on the dust distribution in SMGs, which remains poorly 
known. If the star forming clouds and heavy extinction are 
patchy and distributed over the galaxy, then older stars will 



eventually diffuse out of the heavily obscured regions (per- 
haps on Milky-Way-type time-scales of ~ lOOMyr). In this 
case, we would be able to observe the bulk of the stellar mass, 
except in the earliest formation stages of the galaxy. On the 
other hand, if the entire galaxy is enshrouded in dust, and the 
older stars never diffuse out to regions of lower extinction, a 
large fraction of the stellar mass could be hidden. While com- 
parisons with dynamical mass estimates of SMGs (see §|5J} 
suggest that a large fraction of the stellar mass is not missing 
from our estimates, the dynamical masses are derived using a 
tracer (gas) which may not have the same spatial distribution 
as the stars and may not be representative of entire galaxies. 

5. DISCUSSION 

In this section we combine the stellar luminosities and typ- 
ical stellar mass we have derived for our SMG sample with 
various mass, AGN, and star formation indicators to extract 
their implications for the nature of SMGs and their future 
evolutionary path. We first perform a small, but independent, 
check of our typical stellar mass estimate to evaluate its plau- 
sibility before moving on to perform comparisons of SMGs to 
other populations of high-z galaxies. We will also discuss the 
origin of the non-stellar near-IR continuum emission. Finally, 
we will use constraints from the fading over time of stellar 
population models to predict the z = descendants of SMGs. 

5.1. Comparison to Dynamical Mass from CO Observations: 

Testing M+ 

The various different systematic uncertainties associated 
with stellar mass determinations, and the different results 
which have been obtained in the literature strongly high- 
light the need for independent checks on stellar mass esti- 
mates. In principle, the best check on any stellar mass esti- 
mate i s to compa r e M+ to the dynamical mass, Md y „- How- 
ever, iGreve et all (120051) warn that the dynamical masses of 
SMGs must be treated with caution, because Md y „ for SMGs 
(and other galaxies) is usually calculated assuming the sys- 
tem is in dynamical equilibrium and that we know the ve- 
locity field and mass distribution, whereas SMGs may be the 
product of a merging gas-rich system (e.g..|Smail et alJl!998l 
2004; Tecza et alJl2004t ISwinbank et al.ll2006t iTacconi et alJ 
2008) which is not in dynamical equilibrium. If the sys- 
tem is not in dynamical equilibrium, the mass distribution as- 
sumed is incorrect, or if the observed CO line is not a good 
tracer of the system dynamics or mass, Mdyn will be differ- 
ent. Moreov e r, high-resol ution millimeter obs e rvatio ns by 
ITacconi et aD (120061 120081) and iBothwell etail (1201 Oh sug- 
gest that the CO-emitting region (R ~ l-2kpc) in SMGs 
is frequently smaller t han the overall stellar mass and star 
formation distribution (|Chapman et aT] |2004; Biggs & Ivison 
120081: ISwinban k et alj|2010l) . so the dynamical masses from 
CO observations are often lower limits. A further compli- 
cation is that the dynamical mass is not exclusively a bary- 
onic mass: since a galaxy's dark matter halo contributes 
to the gravitational potential and thus to the velocity field, 
the galaxy's mass derived from the velocity field will con- 
tain some contribution from dark matter. The fractional 
contribution of dark matter to the available CO dynamical 
mass estimates is currently unknown; however, the likely 
descendants of SMGs in the local universe (i.e., galaxies 
of equal or greater stellar mass) have been observed to be 
bary on-dominat ed on scales similar to the CO-emitting region 
in SMGs (e.g.. iGerhard et alJ l200lt ICappellari et all 120061 
iKassin et a l. 2006; Willi ams et alj|2009b . 
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In spite of the uncertainties, we use the 16 SMGs with CO 
observations in our sample to carry out the best test available 
of the stellar masses we have derived. We list the molec- 
ular gas results and dynamical mass estimates for the 16 
SMGs which have published CO observations in Table|4] not- 
ing that only t he CO -detected SMGs have estimates of Mdyn- 
iTacconi et al.l (12008b have carried out a comparison of gas, 
stellar, and dynamical mass for four of the SMGs in our sam- 
ple already, in order to constrain the IMF and CO-to-H2 con- 
version factor for high-z galaxies, and found that the stellar 
masses of two of the SMGs exceeded the dynamical mass. 
However, we are more concerned here in determining if the 
typical stellar mass of SMGs is consistent with the typical dy- 
namical mass, since our stellar masses for individual SMGs 
are highly uncertain. The median dynamical mass for SMGs 
in our sample with CO observations is (Md yn ) = 2.3 x lO n M 
and thus consistent with the sum of the median stellar mass 
(M* = 7.2 x 10 10 M©) and the median molecular gas mass 
[M(H 2 ) = 3.0 x 10 10 M Q ] of the SMGs with CO observations 
(which is identical to the median of th e combined samples of 
iGreve et alj2005t ITacconi et a l. 2006, 2008). Thus, assuming 
the median Mdyn an d MQHz) are representative of the whole 
sample, we conclude that our median stellar mass estimate is 
generally consistent with the published dynamical masses of 
SMGs. 



5.2. The Origin of the Near-IR Continuum Excess in SMGs 

Although the majority of our SMG sample does not con- 
tain a dominant near-IR PL component, we nevertheless 
wish to understand the origin of the near-IR excess ob- 
served in a significant fraction of the galaxies in our sam- 
ple. To test the origin we need to distinguish the near- 
IR emission from the galaxies' nuclei from that which fol- 
lows the obscure d star formation (which may be more dis- 
tributed; see. e.gJChapman et al.l2004tlBiggs & Ivisonl2008t 
Menendez-Delmestr e et al.ll2009l) . Yet, because high-z SMGs 
are unresolved in our IRAC images, we cannot directly con- 
nect the PL emission to the galaxies' nuclei or star formation 
regions; we must find some other means of discriminating be- 
tween the central SMBH and the host galaxy as the source of 
the near-IR excess. One way to do so is to examine other AGN 
indicators in the SMGs: if the origin of the excess is AGN, we 
would expect that the presence and/or luminosity of a signifi- 
cant PL component correlates with other AGN signatures, un- 
less the AGN is highly obscured and Compton-thick (which 
does not appear to be the case for the majority of SMGs; see 
I Alexander e taT1l2005l) . 

The hard X-ray luminosity of a galaxy is regarded as one 
of the best tracers of SMBH activity for Compton-thin nu- 
clei, since the hard X-rays are less susceptible to the heavy 
obscuration thought to enshroud a central SMBH. More- 
over, numerous authors have shown that the hard X-ray lu- 
minosity of QSOs and Seyfert 1 galax ies is correlated with 
the nuclear near-IR luminosity (e.g.. ICarleton et al.l 119871; 
lAlonso-Herrero et alJll997t iMushotzkv et al.ll2008l) . who fre- 
quently conclude that the non-stellar near-IR (e.g., K, L, and 
M bands) emission originates from hot dust heated by the 
AGN. 23 SMGs in our sampl e have rest-frame hard X-ray 
data available in the literature (Alexander et al. 2005, 2008a; 
iLe Floc'h et al. 1 120071) . For these galaxies we can most di- 
rectly test if the near-IR excess is associated with an AGN, 
by checking if the hard X-ray luminosity correlates with the 
luminosity in the near-IR PL component. Although the most 



45 



cn 



44 



43 



42 



41 



CP 

o 




40 



42.5 



Moroston (2005) Models 
CSF SFH 



43.0 



43.5 44.0 
loq L x (erqs/s) 



44.5 



45.0 



Figure 12. Rest-frame A"-band luminosity in the PL component as a function 
of extinction-corrected X-ray luminosity for the SMGs which have X-ray ob- 
servations in the literature. The dotted line represents Ljc(PL Component) = 
Lx . The shaded area shows the region in which local AGN fall. With the ex- 
ception of 2 SMGs which may be Compton-thick AGN, and 3 SMGs which 
have no near-IR PL component, Lx is clearly correlated with the A"-band lu- 
minosity of the non-stellar continuum excess; the correlation suggests that 
AGN are the origin of the near-IR continuum excess in SMGs. 



significant correlations in local AGN samples between near- 
IR and X-ray emission occur in L and M bands due to stel- 
lar contamination in bluer bands, our IRAC data cannot con- 
strain bands redward of rest-frame K for most of our sam- 
ple. Thus, to test for a correlation, we have calculated the 
rest-frame Zf-band luminosity contained within the PL com- 
ponent for each of the 23 SMGs with X-ray data (in add ition 
to computing the same for rest-frame //-band in jj4.11 >. In 
Figure [12] we show the absorption-corrected X-ray luminos- 
ity (Lx) versus the rest-frame /if -band luminosity (Lk) of the 
PL component for the SMGs which have X-ray observations. 
Excepting 2 galaxies with low Lx and high Lk in the PL com- 
ponent, which may host Compton-thick AGN, and 3 galaxies 
which have no near-IR PL component and thus likely lack a 
hot dust component (similar to the local ULIRGs NGC 6240 
and NGC 4945), Lk of the PL component is well-correlated 
with Lx- The correlation remains apparent even when we 
restrict the redshift range of the X-ray observed sample to 
2.0 < z < 2.6, indicating that the correlation is not merely 
an artifact of galaxy distance. The median ratio of the Zf-band 
PL luminosity and Lx for all of the SMGs with X-ray data, 
L K (PL) /L x = 0.9 ± 0.4 with RMS scatter 1 .9, is not dissimilar 
to the Lk/Lx rati o of the average radio-quiet QSO SED from 
lElvis et all (fl994l) (L K /L X ~ 0.4), and the range spans those 
observed for local PG QSO s, Seyfert Is, and Seyfert 2s (e.g., 
lAlonso-Herrero et alj|1997l) . 

The observed L x -L x correlation for our SMGs with X-ray 
data and the similarity of their L x /Lx ratios to local obscured 
and unobscured AGN indicate a probable AGN origin of the 
near-IR excess for the galaxies in which the excess is present 
in the Zf-band. To lend credibility to this conclusion we also 
examine the spectroscopic classifications of the SMGs which 
have significant near-IR PL contributions. When we search 
through the rest-frame UV, Ha, mid-IR, and X-ray spectral 
classifications available for the SMGs in our full sample, we 
find that, of the 25 SMGs which have fractions of rest-frame 
//-band PL-component luminosity greater than 0.2, 20 have 
been classified as AGN. The remaining 5 have starburst classi- 
fications, but the only available spectroscopic data is from the 
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rest-frame UV, which is very sensitive to dust obscuration and 
thus is not a reliable indicator of AGN activitvPl. Thus, the 
spectroscopic classifications of the SMGs in our sample with 
significant PL components are generally consistent with an 
AGN origin for the near-IR continuum excess in those galax- 
ies. However, for the SMGs which are at z < 1.8 and show 
a near-IR excess only in one appropriate photometric band, 
we cannot rule out a contribution from a 3.3 /jm PAH emis- 
sion line without spectroscopic data. We note that SMGs in 
our sample without a near-IR excess do not universally lack 
an AGN spectral identification: in fact, ~ 60% of the galaxies 
without a near-IR excess have been identified as AGN through 
spectroscopy at some wavelength. However, as we mentioned 
in the previous paragraph, some local AGN lack hot dust com- 
ponents; thus, the presence of AGN in SMGs without a near- 
IR excess does not render invalid our conclusion for SMGs 
with significant excess near-IR emission. 

Sinc e such a large fraction of our SMG sample (~ 70%; see 
33. 1 .2b display this near-IR excess which seems to be asso- 
ciated with AGN (at least in the galaxies with the strongest 
PL contribution), it is interesting to determine if the fre- 
quency of near-IR continuum excess observed in our sample 
is similar to that of independent SMG samples. One com- 
pletely independent SMG sample we can compare to is that 
catalogued in the Extended- Chandra Deep Field-South (E- 
CDFS) bv lWeiB et all d2009l) . To reduce the effect of spec- 
troscopic bias from our comparison, we examine the number 
of SMGs in each sample displaying a near-IR excess relative 
to the full number of SMGs in the parent sample, including 
those whose ra d io or mid-IR counterparts were not identified. 
IWardlow et al.l (1201 ll) find that 13 out of the parent sample 
of 126 SMGs in the E-CDFS show a non-stellar excess at 
observed-frame 8.0/im (~ 10%), whereas for our sample 44 
galaxies out of the parent sample of 150 SMGs in C05 dis- 
play non-stellar excess in their IRAC photometry (~ 30%). 
The large difference in frequency of near-IR excess suggests 
that some bias towards hot dust emitters and AGN exists in the 
C05 SMG sample. Some bias toward AGN in the C05 sample 
could have arisen through the spectroscopic selection, in that 
SMGs with multiple emission lines in their spectra were more 
likely to be identified than, for example, galaxies with one (or 
no) emission line. However, future submillimeter surveys in 
additional sky fields are required to confirm that the C05 SMG 
sample suffers from an AGN selection bias. 

5.3. Comparison of SMGs to Other High-z Galaxy 
Populations 

In Paper I we compared the rest-frame near and mid-IR 
colors of SMGs to the samples of stellar-dominated, rest- 
frame-UV-selected, z ~ 3 Lym an-break gala x ies (L BGs) and 
Z ~ 2 BX/BM galaxies from iReddy et all (120061) and the 
radio-loud, obs cured AGN-domina ted high-z radio galaxies 
(HzRGs) from lSevmour et all (120071) . all of which have spec- 
troscopic redshifts. Our comparisons indicated the SMGs are 
brighter and redder than UV-selected galaxies at high-z, sug- 
gesting that they have higher mass, extinction, and possibly 
hot dust contribution. SMGs are similarly bright and red as 
HzRGs, suggestive of similar mass and hot dust contribu- 
tion. Now we have stellar luminosities and a median stellar 

14 As an illustration of the capability of the rest-frame UV spectra to be 
misleading, we note that the two SMGs which do not follow the L\ — Lk 
relation and could host Compton-thick AGN have starburst classifications in 
the rest-frame UV, but are continuum-dominated in their mid-IR spectra from 
IRS. 



mass for the SMG sample and thus we can compare the stellar 
components of the high-z populations directly. As in Paper I, 
we prefer to compare only to galaxy samples with spectro- 
scopic redshifts to eliminate uncertainty associated with red- 
shift; we thus use the same comparison samples of LBGs, 
BX/BM galaxies, and HzRGs as in Paper I. For the present 
study we now also include the z ~ 2 ^-selected galaxies (dis- 
tant red galaxies, or DRGs) wi th spectroscopic redshifts from 
iKrieketalJ (120061 12007, 2008 ), to broaden our range of high-z 
comparison galaxy populations. We note that each sample of 
galaxies representing each population have been selected by 
different criteria and at different wavelengths, and thus each 
suffers from different selection effects which are beyond the 
scope of this paper to analyze. As in Paper I, our compari- 
son is not intended to be comprehensive or representative; we 
merely wish to obtain an idea of the relationships between 
well-studied galaxy samples. 

Looking at the median rest-frame //-band absolute magni- 
tudes for the various types of high-z galaxies, we find that 
the typical DRG, SMG, and HzRG are ~ 2 or more magni- 
tudes brighter than the typical LBG or BX/BM galaxy. The 
reddening-corrected median Mh values for DRGs, SMGs, 
and HzRGs are -26.2 ±0.1 (A. Muzzin, private communi- 
cation), -26.0 ±0.2, and -25.7 (ISeymour et all 120071) . re- 
spectively, wh ile the median Mh for the LBGs and BX/BM 
galaxies in the Reddv et"all d2006l) sample are -23.8 ±0.2 and 
-23.6±0.1, respectively. The ~ 2-2.5 mag difference be- 
tween the UV-selected galaxies and the DRGs, SMGs, and 
HzRGs is statistically significant in all cases where a boot- 
strapped error in the median is available, and corresponds to 
a factor of ~ 6 — 10 in luminosity. The difference in near- 
IR stellar luminosity probably reflects a lower stellar mass in 
LBGs and BX/BM galaxies, for three reasons: (1) the hot dust 
contribution to Mh has been removed from the SMGs and 
HzRGs; (2) the UV-selected galaxies are less obscured, so ex- 
tinction is unlikely to be a factor in the lower Mh', and (3) 
because we receive UV flux from young stars in UV-selected 
galaxies, it is unlikely that the dominant stellar populations in 
LBGs and BX/BM galaxies are older (and thus fainter) than in 
the DRGs, SMGs, and HzRGs. The smaller differences in M H 
among DRGs, SMGs, and HzRGs, however, are less clear in 
origin but are usually not statistically significant. Hence, we 
conclude that DRGs, SMGs, and HzRGs have essentially sim- 
ilar stellar near-IR luminosities, and possibly stellar masses. 

Our next step should be to compare the stellar masses of the 
SMGs, LBGs, DRGs, HzRGs, and BX/BM galaxies. Due to 
the various systematic uncertainties associated with esti- 
mates made through fitting stellar population synthesis mod- 
els, our preference is to compare galaxy samples whose mass 
estimates have been derived in a consistent manner. However, 
no such concordance exists in the literature at this moment 
between the stellar mass-estimating methods of the high-z 
galaxy types selected through different criteria. 

In the absence of consistently-derived stellar mass estimates 
in the literature for different high-z populations, we have cho- 
sen to apply our stellar population model plus power law SED 
fitting techn i que to the sample of UV-selected galaxies in 
Redd y et al.l (120061) and derive stellar masses for that sample 
in the same manner as we have done for the SMGs in our sam- 
pl e. We have chosen to use the LBGs and BX/BM galaxies 
in lReddv et alj ([2006) because complete optical, near-IR, and 
IRAC photometry are readily available for a large number of 
galaxies, along with spectroscopic redshifts. The median stel- 
lar masses we calculate from our method are (M + ) = (7.7 ± 
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0.2) x 10 9 M Q for LBGs and (M*) = (5.9 ±0.1) x 1O 9 M 
for the combined BX/BM sample, whi ch are rough l y con- 
sistent with the stel lar mass estima tes of iReddv et al.l (2006) 
when scaled to the iKroupal (|2001) IMF. When we compare 
the median stellar mass of SMGs, (M*) = 7.2 x 1O 1O M , to 
the self-consistently calculated median stellar masses of the 
UV-selected galaxies, we unambiguously determine that the 
UV-selected galaxies are approximately an order of magni- 
tude less massive than SMGs. Such a disparity in mass does 
not, however, rule out z ~ 2 SMGs as direct evolutionary de- 
scendants of UV-selected galaxies at z ~ 3, because, assuming 
LBGs contain the necessary quantity of molecular gas, the 
time to form the stellar mass of a typical z ~ 2 SMG from 
the typical mass of a typical LBG at the typical LBG star 
formation rate (~ 5OM vr" 1 : iKornei et aill2010h is ~ 1 Gyr, 
roughly the separation in cosmic time between z ~ 2 and 
z~3. 

5.4. The Future Evolution of SMGs 

We have now established the typical stellar content of 
SMGs at z ~ 2; yet we still wish to understand what form 
the descendants of SMGs will have in the local Universe. To 
explore the future evolution of SMGs after the termination of 
the starburst we can employ simple stellar evolution models 
based on these star formation time-scales and use the models 
to predict the stellar luminosity distribution of the descendants 
of SMGs from their Mh distribution at z ~ 2. For the purposes 
of our discussion here, since we could derive no constr aints 
on star formation history from our SED fits (see 33.1. lb . we 
conservatively assume a simple scenario in which we observe 
the z ~ 2 SMGs in the middle of last major starburst event in 
their history; after the burst ends, they will not undergo subse- 
quent mergers or accretion and the stars will evolve passively. 
Our scenario thus provides a lower limit on the likely z = 
luminosity of SMG descendants. 

We first construct the high-z Mh distribution for our sample 
of z > 1.5 SMGs using the power-law continuum-subtracted 
(stellar) values of Mh which have been individually de- 
reddened and averaged between the IB and CSF star forma- 
tion histories. The final distribution, shown in Figure [TH as 
the solid histogram, has median (Mh) =—26.0 ±0.2. We note 
that we do not exclude from the stellar Mh distribution the 
SMGs which have significant (> 50%) PL components, al- 
though their stellar Mh values probably have large uncertain- 
ties. Then, we evolve our sample of z > 1.5 SMGs to the 
prese nt day using the M05 stellar population synthesis models 
with a Kroupa (2001) IMF. We consider as possible scenarios 
that SMGs are starbursts of constant SFR of duration 50, 100, 
200, and 400 Myr, covering the range of estimates of the dura- 
tion of the starburst in the SMG phase from CO observation s 
(e.g. JNeri et al.ll2003HGreve et al.ll200l iTacconi et al.ll2006h . 
We assume that we observe each SMG halfway though its star 
formation burst, and then compute the fading that occurs be- 
tween ages 25, 50, 100, and 200 Myr and 10 Gyr (roughly the 
difference in cosmic time between z ~ 2 and z = 0). The faded 
luminosity distributions for each burst duration are shown as 
the dotted histogram in Figures [13] and [14] and predict me- 
dian stellar absolute H magnitudes for the SMGs at z = of 
Mh = -22.1,-22.6,-23.1, and -23.6, respectively, for burst 
durations of 50, 100, 200, and 400 Myr. We observe that the 
predicted me dian values of Mh are all fainter than L* at z = 
(M* H = -23.7; ICole et al.ll2001l) ; however, recall that these pre- 
dictions are conservative lower limits. Subsequent mergers or 



star formation can increase the mass and luminosity of SMG 
descendants by varying factors, and may reasonably be ex- 
pected. 

We can put these aged and faded SMGs into context by 
comparing them to a sample of massive elliptical galaxies in 
the local Universe. We use as our local massive elliptical 
galaxy sample the ellipticals in the Coma cluster since they 
have a large, homogeneous, near-IR photometry data set and 
optical morphology data set. We construct the Mh distribu- 
tion for massive ellipticals in Coma using //-band 2MASS 
photometry and optical morphologies compiled by R. Smith 
and J. Lucey (2009, private communication). To convert the 
observed H magnitude s to Mh, we use a d istance modulus 
m-M = 35.0mag (e.g. jLiu & Grahamll2001l) . We include the 
Mh distribution of Coma ellipticals in Figures [13] and [14] as 
the dashed histogram. Somewhat surprisingly, the SMG dis- 
tributions faded for 50, 100, and 200 Myr bursts all reason- 
ably match the Coma Mh distribution within ~ 0.5 mag: two- 
sided Kolmogorov-Smirnov (K-S) tests comparing to Coma 
distribution to each of the faded SMG distributions produce 
probabilities of being drawn from the same parent distribu- 
tion of 1.3%, 35%, and 13%, respectively, so we cannot dis- 
tinguish any of these faded distributions from the Coma el- 
lipticals. (We note that the 400 Myr burst distribution, when 
compared to the Coma distribution in the K-S test, produces a 
probability of 0.07%.) As shown in Figure[l3] the SMG distri- 
bution faded according to the 100 Myr burst model appears to 
best match the distribution of the Coma cluster ellipticals, al- 
though the differences between the 50, 1 00, and 200 Myr burst 
models applied to the SMG distribution are probably not sig- 
nificant. Our comparison of the Mh distribution of the ellipti- 
cal galaxies in the Coma cluster to faded Mh distributions for 
Z > 1 .5 SMGs thus suggests that, for reasonable starburst du- 
rations (50-200 Myr), and under the conservative assumption 
of no significant growth in stellar mass after the SMG phase 
terminates, the z = SMG descendants may be massive ellip- 
tical galaxies with a luminosity (and therefore mass) distri- 
bution similar to that of massive ellipticals in the Coma clus- 
ter (the majority of which fall in the range Q.2L* — 2L*). We 
note that this range is somewhat lower than t he luminosity of 
presen t-day SMG descendants estimated by Swinbank et al. 
(2006, > 3L*) by comparing the average velocity dispersion 
of 5 SMGs to the present-day Faber-Jackson relation. In light 
of the uncertainty in both analyses (arising from star forma- 
tion history assumptions in this paper and the assumptions of 
virial equilibrium and no evolution in velocity dispersion by 
Swinbank et al.), the two results may be reconcilable. How- 
ever, our sample of S MGs is ~ 10 times larger than that of 
ISwinbank et alj d2006l) . and therefore is likely more represen- 
tative of the range of mass and luminosity among SMGs. 

If burst durations of 50-200 Myr are indeed correct for 
SMGs, then it is interesting to determine from them what frac- 
tion of today's massive ellipticals underwent a SMG phase. 
We can use the starburst durations which predict the SMG 
fading that best matches the Coma ellipticals to estimate 
the space density of SMG descendants of similar luminos- 
ity/mass to Coma ellipticals, and see how that compares to 
the space density of similar ellipticals from local galaxy lumi- 
nosity functions. We calculate the space density of SMG de- 
scendants assuming each galaxy goes through only one SMG 
phase using the simple relation 

<t>D = PSMG (fobs/ fburst)j (2) 

where (j) D is the comoving space density of the descendants 
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Figure 13. Comparison of the absolute H magnitude (Mm ) distribution of the elliptical galaxies in the Coma cluster to that of the stellar component of our 
spectroscopically-identified SMGs at z > 1.5. The distributions match well if we passively evolve the stars in SMGs to an age of lOGyr, modeling the SMGs as 
100-Myr-long starburst events which, at z ~ 2, we observe halfway through at an age of 50 Myr. The match between the faded SMG distribution and the present- 
day Coma distribution suggests that, for the predicted time-scale of the starburst phase in SMGs, SMGs will evolve into galaxies of similar stellar luminosity and 
mass as the massive ellipticals in the Coma cluster if they evolve passively upon termination of the starburst. 



of SMGs, psmg is the observed space density of SMGs, t a b s 
is the epoch over which we observe SMGs, and fb urs t is the 
duration of the starburst. We determine f t, s from the redshift 
distribution for our SMG sample (updated fro m C05 with new 
results from Menendez-Delme stre et al.l2 009) after correcting 
fo r spectroscopic and radi o detection incompleteness; similar 
to Swinba nk et al.l (|2006), we use the ±lrr interval around 
the median of the redshift distribution to define t b s — 2.3 Gyr. 
We estimate the observed space density of SMGs from the 
integral submm source counts for S&so > 4mJy found by 
ICoppin et ail d2006l) . 844 ± 1 17 deg" 2 , and the comoving vol- 
ume between z = 1.46-3.06. In this way, the 50, 100, and 
200 Myr starburst durations imply that the space density of 
SMG descendants at z = will be <f> D ~ 2 , 1 , 0.5 x 1 0" 3 Mpc" 3 . 
When we integrate the ^T- band luminosity fu nction of early- 
type galaxies at z = of ICroton et al.l d2005l) . we find that 
we can match the space density of SMG descendants of 
4>d ~ 0.5 x 10~ 3 Mpc~ 3 by integrating the early-type luminos- 
ity function down to 1.3L*. For the lOOMyr-duration burst, 
we can match the predicted space density of 1 x 10~ 3 Mpc~ 3 
if we integrate the early-type luminosity function down to 
0.8L*. For the case of the 50 Myr-duration burst, we can 
match the predicted space density of 2 x 10~ 3 Mpc~ 3 by in- 
tegrating the early-type luminosity function to 0.4L*. Thus, 
our calculations suggest that for starburst durations of 50- 
200 Myr, SMGs could account for the formation of the entire 
population of elliptical galaxies with luminosities greater than 
-0.4L*-1.3L*. 

Given the systematic uncertainty associated with the form 
of the stellar IMF for high- z gal axies, our median stellar mass 
for z > 1 .5 SMGs from §14.31 is essentially consistent with 
m*, and is consistent with scenarios which predict that SMGs 
evolve into massive elliptical galaxies. Moreover, the bary- 
onic masses (gas plus stars) of high-z SMGs, which serve as 



an upper limit on the final stellar mass of SMGs assuming 
no future accretion or gas loss, also point toward a future for 
SMGs as ~ m* galaxies. Assuming t he m edian molecular gas 
mass of CO-observed SMGs (from ^5. Il l is typical of our en- 
tire sample in an average sense, then the typical baryonic mass 
(gas plus stars) for our SMGs is (M barvon ) ~1.0x 10 n M Q . 
This estimate excludes the contribution from neutral hydro- 
gen; if we reasonably estimate that the contribution of neu- 
tral gas is similar to that in the Milky Way and thus equal to 
the mass of molecular gas, then the typical baryonic mass of 
SMGs would be (M harvon ) ~1.3x 1O H M . Both estimates 
of baryonic mass closely resemble the stellar mass of a L* 
galaxy in the local universe (1.4 x 10"M ? ; ICole et al.ll2001l) . 
On the other hand, feedback from star formation and/or AGN 
may result in outflows of gas from the galaxy so that only a 
fraction of the gas in SMGs at z ~ 2 is eventually converted to 
stars (e.g., Begel mar][l98l ISilErT997t iDiMatteo et aHl200l 
Murray et al.l 120101) . Nevertheless, the baryonic mass is a 
strong indicator that the end product of the starburst phase 
in a typical SMG will be a galaxy of mass of order m* at 
z = 0, within a factor of a few. Even considering the uncer- 
tainties, the baryonic mass is too small for typical SMGs to 
represent the formation phase of the very luminous, rare, cD- 
type galaxies observed in galaxy clusters; the space density of 
SMGs is also far too high to expect typical SMGs to all be so 
massive. As shown in Figure [10] a range of approximately an 
order of magnitude in stellar luminosity and mass does exist 
among z ~ 2 SMGs, including a high-mass tail, so some SMG 
descendants wil l evolve into ma s sive el lipticals of 3L* and up 
as predicted by Swinban k et al.l (120061) . Subsequent mergers 
or gas accretion will also produce more massive descendants. 
However, in the absence of significant mass accretion, the av- 
erage SMG seems destined to be an ~ L*-type galaxy. 

6. SUMMARY 
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Figure 14. Comparison of the absolute H magnitude (Mm ) distribution of the elliptical galaxies in the Coma cluster to that of the stellar component of our 
spectroscopically-identified SMGs at z > 1.5, assuming starburst lifetimes for the SMGs other than 100 Myr. The fading predicted by modeling SMGs as 
starbursts of duration 50 and 200 Myr cause the Mh distribution of SMGs to match the Coma cluster similarly as the 100 Myr burst; however, a 400 Myr burst 
duration will produce SMG descendants which are more luminous than is typical of the massive ellipticals in Coma. 



In this paper we have fit the rest-frame UV through near- 
IR SEDs o f the radio-detec t ed SM Gs with spectroscopic red- 
shifts from Chapma n et al.l d2005) to estimate stellar masses. 
For the first time we remove the excess continuum contribu- 
tion to the galaxies' near-IR luminosity to ensure that we de- 
rive a mass for the stellar component only. Observing that the 
K-band luminosity of this near-IR excess is well-correlated 
with hard X-ray luminosity of the galaxy, we conclude that 
the near-IR continuum excess is caused by the galaxies' AGN; 
however, for only 1 1 % of our SMGs is the AGN contribution 
stronger than the stellar contribution in the rest-frame near- 
IR. In trying to fit synthetic stellar populations to our contin- 
uum excess-subtracted SEDs, we show that our photometric 
data spanning a broad range of wavelengths are still not suffi- 
cient to place useful constraints on the star formation history 
of the galaxies and, thus, the age of the dominant stellar pop- 



ulations, since our galaxies are heavily reddened and our data 
are not well-enough sampled near the strong age indicators of 
the Balmer and 4000 A breaks. 

The median stellar mass of the galaxies in the sample is 
better constrained, however, despite systematic uncertainties. 
By taking a simple, physically motivated approach using a 
constant mass-to-light ratio, we find that SMGs have a mod- 
est median stellar mass of 7 x 1O 1O M0, which is lower than 
most previous estimates of stellar mass in SMGs. Our new, 
lower stellar mass estimate for typical SMGs is consistent 
with the median dynamical mass of SMGs with CO obser- 
vations. When we compare our typical SMG stellar mass to 
those of other high-z galaxy populations, the SMGs appear 
to be more massive by a factor of 10 than high-redshift UV- 
selected galaxies. 

We then explore the possible evolution of SMGs after they 
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run out of molecular fuel for star formation, based on es- 
timates in the literature for the gas consumption/starburst 
timescale in SMGs. We find that if we passively evolve SMGs 
to z = 0, the range of expected starburst time scales causes the 
rest-frame H absolute magnitude distribution of high-z SMGs 
to fade so that it nearly matches the M H distribution of mas- 
sive elliptical galaxies in the Coma galaxy cluster. The match 
to the Coma cluster distribution suggests that the descendants 
of typical SMGs will be <~ m* galaxies similar to Coma el- 
lipticals in the local universe, and is further supported by the 
typical baryonic mass of our SMG sample (~ 10 h Mq) which 
places an upper limit on the future stellar mass of SMGs as- 
suming no further cold gas accretion. Thus, our results con- 
firm the picture of high-z SMGs as highly luminous, massive 
galaxies simultaneously experiencing strong, obscured star- 
burst activity and AGN activity fueled by massive molecular 
gas reservoirs, which will evolve into ~ L*-type galaxies fol- 
lowing the exhaustion of their molecular gas. 
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Table 1 

Photometric Data Used in SED Fits of C05 SMGs 



Field Name 


Optical Bands 


Near-IR Bands 


IRAC Bands 


{# Optical Bands)" 


Galactic E(B-V) b 








(fim) 


(per galaxy) 


(mag) 


CFRS-03h 


U,B,V,R,I 


J,K 


3.6,4.5,5.8,8.0 


6 


0.071 


Lockman 


B,R,I 


K 


3.6,4.5,5.8,8.0 


3 


0.007 


GOODS-N 


U,B,V,R,I,z' 


J,K 


3.6,4.5,5.8,8.0 


8 


0.011 


SSA-13 


B,R,I,z 


J,K 


3.6,4.5,5.8,8.0 


6 


0.014 


CFRS-14h 


B,R,I 


J,K 


3.6,4.5,5.8,8.0 


5 


0.011 


ELAIS-N2 


B,V,R,I 


K 


3.6,4.5,5.8,8.0 


5 


0.006 


SSA-22 


B,R,I 


J,K 


3.6,4.5,5.8,8.0 


5 


0.061 



The median number of optical and near-IR photometric bands per SMG actually used in the SED fits. 
Estimated line-of-sight reddening due to the Milky Way from Schlegel et al. ( 1998). 



Table 2 

Results of SED Fitting and Stellar Masses for C05 Spectroscopic SMG Sample, M05 

Models 



Chapman et al. (2005) ID 


Zspcc 


M H (obsf 


A///(PL-sub) b 


a c 


PL Fraction (IB) d 
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2.689 


> -26.29 


> -26.29 


2.0 


q 0+0.6 


0.0+ 05 


1.8 ±0.3 


< 10.93 


SMMJ105219.15+571858.4 


2.372 


-26.28 ±0.01 


-26.28 ±0.01 


2.0 


0.0+ 04 


0.0+ 04 


1.5±0.0 


10.93 ±0.05 
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3.0 
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-25.01 ±0.04 


2.0 


°- 8+ 02 
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10.42 ±0.05 
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2.509 


-26.95 ±0.01 


-26.95 ±0.01 


2.0 


0.Q +M 
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11.20 ±0.05 


SMMJ123616.15+621513.7 


2.578 


-26.91 ±0.01 


-25.84 ±0.73 


2.0 


9+° 1 


6+° 2 


1.4 ± 1.4 


10.85 ±0.27 


SMMJ123618.33+621550.5 


2.000 


-25.85 ±0.07 


-25.75 ±0.07 


3.0 


5+ 02 




1.3±0.4 


10.71 ±0.05 


SMM J 1 2362 1 .27+62 1 708.4 


1.988 


-25.77 ±0.03 


-25.64 ±0.08 


3.0 


04+0.1 




1.0 ±0.4 


10.67 ±0.05 


SMMJ123622.65+621629.7 


2.466 
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3.0 


2+^ 
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6+ 8 - 1 
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-26.29 ± 0.07 


-26.02 ±0.06 
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10.82 ±0.04 
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-27.90 ±0.14 
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3.0 


2.2 ±0.5 


11.29 ±0.10 
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-23.54 ±0.04 


-23.52 ±0.05 


2.0 
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9.83 ±0.06 
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-26.93 ±0.01 


-26.93 ±0.01 


2.0 
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2.565 


-28.32 ±0.10 


-27.34±0.18 


3.0 


0.4 ±0.4 


11.36±0.10 


SMM J 1 3 1225 .20+424344.5 


1.038 


-25.73 ± 0.04 


-25.60 ±0.04 


3.0 


1.4±0.1 


10.66 ±0.04 


SMMJ131225.73+423941.4 


1.554 


-25.25 ±0.10 


-25.25 ±0.10 


2.0 


0.0+ 05 


0.0|; 


1.0 ±0.4 


10.52 ±0.05 


SMM J 1 3 1 228. 30+424454. 8 


2.931 


-26.65 ± 0.02 


-26. 14 ±0.10 


3.0 


6+ 01 
u.o 2 

l-O-o.i 


7 +0.f 

u -'-0.2 

l-O-o.i 


1.2±0.0 


10.87 ±0.05 


SMMJ131232.31+423949.5 


2.320 


-27.68 ±0.01 


-24.31 ±0.03 


2.0 


4.8 ±0.0 


10.14±0.05 


SMMJ131239.14+424155.7 


2.242 


-25.44 ±0.06 


-25.41 ±0.09 


2.0 


0.0+ 03 


0.2+ 05 

9+^ : 1 
-0.1 


1.3±0.4 


10.58 ±0.05 


SMM J141741. 81+522823.0 


1.150 


-28.53 ±0.06 


-28. 12 ±0.09 


3.0 


9+° 1 


2.8 ±0.4 


11.67 ±0.07 
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Table 2 — Continued 



i^napman ei ax. ) id 


Zspec 


A4 tnUcY<l 


A/f /PI cuMb 


c\ 


PL Fraction (IB) d 


r l Fraction ^Ljr ) 


Ay 
\ LLLa &/ 


10gM* s 
(Mr,) 


SMMJ141742.04+523025.7 


0.661 


-25.66 ±0.15 


-25.49 ±0.15 


2.0 


0.SJJ" 
0.0 +ii:i 


0.8^-J 


0.6 ±0.3 


10.62 ±0.09 
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2.128 


> -23.65 


> -23.65 


2.0 


0.0 +ii:4 




< 9.88 


SMM J 14 1 800.40+522820.3 
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-26.67 ±0.09 


-26.56 ±0.10 


3.0 


5 +01 
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1.8±0.3 


11.05 ±0.07 


SMMJ141802.87+523011.1 


2.127 


-24.30 ±0.01 


-24.27 ±0.02 


3.0 


3+ 05 


0.0+ 03 


0.2 ±0.2 


10. 12 ±0.05 
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2.712 
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2.0 
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11.53 ±0.04 
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2.283 


-26.19 ±0.07 


-25.84 ±0.09 


2.0 


6 +0 2 
9 +0 ' 


1.0 ±0.4 


10.75 ±0.05 


SMMJ163639.01+405635.9 


1.495 


-25.74 ±0.06 


-25.43 ±0.11 


3.0 


1.5±0.6 


10.59 ±0.07 


SMM J 1 63650.43+405734.5 


2.378 


-27.07 ±0.08 


-26.39 ±0.28 


2.0 




if!-? 

6 +0T 
uo -0.4 


0.6 ±0.3 


10.99±0.13 


SMM J163658. 19+410523.8 


2.454 


-26.61 ±0.08 


-26.48 ± 0.09 


2.0 




3 +0 - 2 

i?A 


1.6±0.4 


11.01 ±0.05 


SMMJ163658.78+405728.1 


1.190 


-26.22 ±0.06 


-26. 17 ±0.08 


3.0 


04 82 


2.1±0.3 


10.89 ±0.06 


SMM J 1 63704. 34+4 10530.3 


0.840 


-23.57±0.12 


-23.46 ±0.23 


2.0 


0.8* 


0.0 +M 


1.3±0.4 


9.81 ±0.11 


SMMJ163706.51+405313.8 


2.374 


-27.08 ±0.14 


-26.76 ±0.21 


3.0 


0.7+ 01 

6+^ 
-8? 
6+< it 


4+o.i 
0J H 


2.1±0.3 


11.13±0.11 


SMM J221733. 12+001 120.2 


0.652 


-25.60 ±0.10 


-25.53 ±0.04 


3.0 


1.8±0.0 


10.63 ±0.05 


SMMJ221733.91+001352.1 


2.865 


-26.85 ± 0.07 


-26.49 ± 0.04 


3.0 


3 + <> :2 


0.9 ±0.3 


11.01 ±0.05 


SMM J221735. 15+001537.2 


3.098 


-25.71 ±0.07 


-25.22±0.18 


3.0 


3+ 02 
-83 


04 +o.2 


0.8 ±0.5 


10.51 ±0.07 


SMM J221735. 84+001558.9 


3.089 


-26.20 ±0.01 


-25.79 ±0.20 


2.0 


3+^ 


1.2±0.3 


10.74±0.10 


SMM J22 1 804.42+002 1 54.4 


2.517 


-26.18 ±0.01 


-25.92 ±0.03 


3.0 


6+^ 


u -'-0.5 


1.2±0.0 


10.79 ±0.05 


SMMJ221806.77+001245.7 


1.910 


-26.17 ±0.09 


-26. 17 ±0.09 


2.0 


0.0+ 01 


0.0+ 01 


2.1 ±0.3 


10.89 ±0.07 



Absolute rest- frame H magnitude derived from observed photometry, averaged between instantaneous burst and eontinous SFH model fits. 
Absolute rest-frame H magnilude when besi-fil PL fraction is subtracted, averaged between instantaneous burst and eontinous SFH model fits. 

Choice of PL index resulting in best fit for galaxy (lowest x^)- in me cases where the best-fit PL fraction is zero, the value listed in this column indicates the index for which the error bar is calculated. 
Best-fit fraction of observed 8.0 /im flux contributed by PL component for instantaneous burst SFH models. 
Best-fit fraction of observed S.O /.iin flux contributed by PL component for continuous SFH models. 

Average of best-fit V-band extinctions from instantaneous burst and continuous SFII lils. The listed error is calculated as half the difference between the two averaged values. 
Stellar mass derived as described in |3.2| The listed error represents half the difference between the instantaneous burst and continuous SFH values. 

SED fits indicate that C05 redshifl of - = 0.276 is incorrect. Photometric redshifl obtained from HYPER-Z suggests alternate redshift of - = 1.81. M* has been calculated using z = 1.81. Due to uncertainty in the redshifl, the results derived for this 
galaxy have not been included in the analysis here. 

i 

SLD I'll miL 1 l.ii Is to converge lor any model, possihh due tu incorrect L'ala\_\ icthhili nr error in couniciparl select inn v, lien nica-.ii nil!: optic, il and IK. -\C plioii unci r>. I lie result lur this enlaxy have not been included in any ol [lie analysis in the le.M. 
^ SMG displays UV/blue excess. 



Known optical QSO. 



Table 3 

Near-IR Continuum Excess-Dominated SMGs" 



Chapman et al. (2005) ID % L H in PL Component (IB) % L H in PL Component (CSF) 



SMM J 1 23600. 1 5+62 1 047.2* 

SMMJ123606.72+621550.7 

SMMJ123616.15+621513.7 

SMMJ123632.61+620800.1 r 

SMMJ131215.27+423900.9 

SMMJ131222.35+423814.1 

SMMJ131232.31+423949.5 

SMM Jl 63650.43+405734.5 



67 
56 
71 
67 
45 
60 
97 
53 



56 
35 
67 
51 
54 
97 
28 



a 
b 
c 



"Continuum-Dominated" is denned as having greater than 509! of the lolal Lft contributed by the PL component. 
Possible Complon-lhickAGN; see Fig. 4 of Alexander el al. (2008b). 
SMG displays UV/blue excess. 
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Table 4 

SMGs with CO Observations 



Chapman et al. (2005) ID z spe c CO Transition M(H 2 ) a M dm Reference 

xl0 10 M Q xIO^Mq 



SMMJ105230.73+572209.5 


2.611 


3-5>2 


< 1.4 




1 


SMMJ105238.30+572435.8 


3.036 


3^2 


< 3.0 




1 


SMMJ123549.44+621536.8 


2.202 


3^2 


3.2 


0.94 


2 


SMMJ123618.33+621550.5 


2.000 


4^3 


7.2 


2.2 


3 


SMMJ123634.51+621241.0 


1.219 


2->l 


6.8 


7.0 


4 


SMMJ123707.21+621408.1 


2.490 


3^2 


2.3 


2.4 


2 


SMMJ12371 1.98+621325.7 


1.992 


4^3 


3.5 




3 


SMMJ 1 3 1 20 1 . 1 7+424208. 1 


3.405 


l->0 


16 


1.2 


1,5 


SMMJ131222.35+423814.1 


2.565 


3^2 


1.0 


0.25 


6 


SMMJ131232.31+423949.5 


2.320 


3^2 


<2.4 




1 


SMMJ163631.47+405546.9 


2.283 


3^2 


<0.9 




1 


SMMJ163639.01+405635.9 


1.495 


2^1 


< 1.8 




1 


SMMJ163650.43+405734.5 


2.380 


3^2 


5.4 


3.4 


2 


SMMJ 16365 8. 19+410523.8 


2.454 


3^2 


4.4 


1.4 


2 


SMMJ163706.51+405313.8 


2.374 


3^2 


2.4 


3.4 


1 


SMMJ221735. 15+001537.2 


3.098 


3^2 


3.0 


2.8 


1 



Note. — References: (1) Greve et al. (2005); (2) Tacconi et al. (2008); (3) Bothwell et al. (2010); (4) Frayer et al. (2008); (5) Hainline et 
al. (2006); (6) Coppin et al. (2008) 

a M(H2) calculated assuming a CO-to-H2 conversion factor a = 0.8 and ^co(4^3)/^co(i^o) = 1-0, ^co<3^2)/^co(i^0) = 1-0, or 

^CO(2^1)/^CO(1^0) = 10- 



